The Astronomy and Astrophysics Review manuscript No. 

(will be inserted by the editor) 



Fundamentals of collisionless shocks for astrophysical 
application, 2. Relativistic shocks 

A M Bykov and R A Treumann* 



Received: date / Accepted: date 

Abstract. In this concise review of the recent developments in relativistic shock theory in 
the Universe we restrict ourselves to shocks that do not exhibit quantum effects. On the other 
hand, emphasis is given to the formation of shocks under both non-magnetised and magne- 
tised conditions. We only briefly discuss particle acceleration in relativistic shocks where 
much of the results are still preliminary. Analytical theory is rather limited in predicting the 
real shock structure. Kinetic instability theory is briefed including its predictions and limita- 
tions. A recent self-similar relativistic shock theory is described which predicts the average 
long-term shock behaviour to be magnetised and to cause reasonable power law distribu- 
tions for energetic particles. The main focus in this review is on numerical experiments on 
highly relativistic shocks in (i) pair and (ii) electron-nucleon plasmas and their limitations. 
These simulations do not validate all predictions of analytic and self-similar theory and so 
far they do not solve the injection problem and the self-modification by self-generated cos- 
mic rays. The main results of the numerical experiments discussed in this review are: (i) 
a confirmation of shock evolution in non-magnetised relativistic plasma in 3D due to ei- 
ther the lepton-Weibel instability (in pair plasmas) or to the ion-Weibel instability; (ii) the 
sensitive dependence of shock formation on upstream magnetisation which causes suppres- 
sion of Weibel modes for large upstream magnetisation ratios a > 10^^; (iii) the sensitive 
dependence of particle dynamics on the upstream magnetic inclination angle 0b„, where par- 
ticles of Obk > 34° cannot escape upstream, leading to the distinction between 'sub-luminal' 
and 'super-luminal' shocks; (iv) particles in ultra-relativistic shocks can hardly overturn the 
shock and escape to upstream; they may oscillate around the shock ramp for a long time, so 
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to speak 'surfing it' and thereby becoming accelerated by a kind of SDA; (v) these particles 
form a power law tail on the downstream distribution; their limitations are pointed out; (vi) 
recently developed methods permit the calculation of the radiation spectra emitted by the 
downstream high-energy particles; (vii) the Weibel-generated downstream magnetic fields 
form large amplitude vortices which could be advected by the downstream flow to large 
distances from the shock and possibly contribute to an extended strong field region; (viii) 
if cosmic rays are included, Bell-like modes can generate upstream magnetic turbulence at 
short and, by diffusive re-coupling, also long wavelengths in nearly parallel magnetic field 
shocks; (ix) advection of such large-amplitude waves should cause periodic reformation of 
the quasi-parallel shock and eject large amplitude magnetic field vortices downstream where 
they contribute to turbulence and to maintaining an extended region of large magnetic fields. 

Keywords CoUisionless shocks ■ relativistic shocks ■ generation of magnetic fields ■ Weibel 
modes ■ Bell modes • Gamma Ray Bursts ■ Pulsar Wind Nebulae termination shocks ■ 
external shocks ■ internal shocks ■ particle acceleration ■ shock radiation ■ downstream 
turbulence 

PACS 95.30.Qd ■ 52.35.Tc ■ 52.72.h-v 
1 Introduction 

Key observations in astrophysics suggest that relativistic collisionless shocks play an im- 
portant if not a central role in the universe. Since many of the astrophysical processes are 
quite violent, i.e. release large amounts of energy on short time scales, it is quite reasonable 
to expect that shocks do frequently form either as explosion driven blast waves, caused by 
ultra-relativistic outflows of matter encountering an obstacle or interacting with other flows, 
or by nonlinear growth and steeping processes of unstably excited waves in those flows; 
and that a substantial fraction of these shocks will be relativistic. These two ways of shock 
formation already lead to the distinction of two types of relativistic shocks: external shocks 
produced in the interaction of flow with an external medium, and internal shocks which 
evolve in the absence of any external obstacle inside the flow. 

Formation of shocks is a well-proven fact in the manifestly collisionless plasma of inter- 
planetary space which has been, is and for long time will be the only place in the Universe 
where large-scale collisionless shocks can be studied in situ. Among others, one important 
lesson that can be learned from their observation is that none of them form on the scale of 
the mean free path A,,,/,, (or 'resistive scale') which in interplanetary space with its character- 
istic dimension ~ 10^ AU to 10^ AU (the size of the heliosphere) is of the order of several 
AU. Without any exception the widths of collisionless shocks are of the order of the ion 
inertial or ion gyro-radius scales, i.e. < XQ^^Xmfp- This implies that even in astrophysical 
systems whose global dimensions L ^ X,„fp grossly exceed their internal mean free paths, 
collisionless shock formation proceeds on micro-scales As that are far below the coUisional 
mean free path with shock transition widths As <S ?i,„fp <C L. 

There are no relativistic shocks in the reach of any man-made spacecraft, however. The 
shocks in the solar system, the heliosphere, and the heliospheric surrounding are without 
any exception non-relativistic. Therefore, all information on relativistic shocks other than 
observational is based on theoretical arguments, ingenious though sometimes questionable 
speculation, numerical simulations or generalisation of our heliospheric knowledge of non- 
relativistic collisionless shocks to the relativistic domain. The most contemporary and com- 
prehensive accounts of the latter have been given in |Balogh and Treumann|p011| l. Such ar- 
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guments rely solely on indirect observations of either high-energy particles (Cosmic Rays), 
intense electromagnetic radiation emitted from certain extended objects, particularly from 
those that are narrow in their transverse dimension, or on Gamma Ray Bursts (GRBs). For 
these reasons nothing about the very structure of collisionless relativistic shocks can be said 
that is comparably precise to our knowledge on non-relativistic shocks. 

From an astrophysical point of view lack of the precise structure of relativistic shocks 
might, in fact, not be as important as it is for the near-Earth non-relativistic shocks. In 
most cases it suffices to assume that certain global conditions are specified at a relativistic 
shock: its speed, Mach number, compression ratio and magnetic-field strength. This is also 
the historical path along which reference to relativistic shocks has in the past been gone. 
With these characteristics in mind that are based on global plausibility considerations, the 
shock serves as an instrument for the generation of the two primary astrophysical effects that 
can be observed by remote sensing: Cosmic Rays and radiation, respectively, energetic (i.e. 
accelerated) particles, and photons. The former are described by a set of measured numbers - 
composition, abundance, anisotropy, and the properties of their energy spectrum like shape, 
slope, breaks, cut-offs etc. -, which contain information about the source of Cosmic Rays, 
their generators and a few properties of the generators. 

Radiation in the electromagnetic spectrum provides information on the spatial location 
of the shocks, their geometrical shape, their relation to the astrophysical object, and on 
processes inside and in the environment of the shock and thus on the mechanism of shock 
production and possibly even the shock structure under various conditions. It is mostly the 
generation of electromagnetic radiation that requires knowledge about the processes that 
generate and maintain shocks. This, however, is a wide field as the generation of radiation 
implies the presence of particles which emit radiation. At a shock, such particles must be 
prepared by the shock to generate the radiation, which is a most complicated and unresolved 
problem. Unlike Cosmic Rays, which reach our solar system, the shock-prepared particles 
cannot be directly observed. 

In order to understand the production of Cosmic Rays, on the other hand, the mechanism 
of shock generation and shock structure are of lesser importance as long as the fraction of 
shock-accelerated energetic particles that are present in the shock environment remains low. 
When this condition is violated, then the Cosmic Rays do also contribute to shock structure 
and formation by "mediating" the shock. 

Because they are very strong radiators (much stronger than non-relativistic shocks) and 
therefore must be efficient particle accelerators, emissions from narrow regions in astro- 
physical space has led to the suspicion that very high Mach number shocks reaching into the 
relativistic domain cause the observed effects. In most of these cases the existence of rel- 
ativistic shocks in astrophysics has been assumed and the conditions from non-relativistic 
shock calculations like, for instance, the shock compression ratio, have been extrapolated 
into the relativistic domain. Though this naturally is a first reasonable approach to an unex- 
plored field, more recent numerical investigation based on particle-in-cell (PIC) simulations 
demonstrated that relativistic shocks exhibit particular properties which are quite different 
from those of ordinary coUisional gasdynamic and also collisionless plasma shocks. These 
properties are related to the self-consistent generation of magnetic fields, production of fast 
particles, which tend to modify the shock, and intense radiation. The present review of some 
of the contemporary achievements extends a previous publication ( [Treumann|2009[ l into the 
relativistic domain. 
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Fig. 1 Sketch of a pulsar-wind shock 
system. The figure is drawn in the 
pulsar system where the supernova- 
remnant (SNR) material streams 
against the ultra-relativistic outflow 
from the pulsar. The interaction be- 
tween the wind and the SNR pro- 
duces a bow shock wave behind 
which a turbulent transition region is 
located. The outer boundary of the 
pulsar wind is a discontinuity, either 
contact or tangential depending on 
the local conditions and on the mag- 
netic field. Inward one expects the 
turbulent pulsar wind transition be- 
tween the discontinuity and the ter- 
mination shock of the wind (after |van"] 
[deTSwaluw et al|2003l . 

2 Observational Evidence for Relativistic Shocks 

Inference on relativistic shocks is obtained either from observations of high-energy radia- 
tion produced by the shock-accelerated particle component or, indirectly, from the detection 
of high-energy Cosmic Rays of non-solar origin. When radiation accumulates around the 
shock, the emitted radiation can be taken as a signature, mapping the shock location and 
for extended shocks also their geometry. Since radiation in the high-energy range comes 
from accelerated particles, all relativistic shocks are in one or another way involved into 
the acceleration of charged particles. It is intuitive that the efficiency of acceleration should 
drastically increase with shock Mach number and thus be reflected in the hardness of the ac- 
celerated particle spectrum and the highest achieved particle energies. Moreover since syn- 
chrotron emission in magnetic fields and bremsstrahlung are closely related to the hardness 
of the particle spectrum, observation of flat radiation spectra and radiation in the Gamma 
and X-ray range can be taken as indication of high relativistic Mach numbers and the pres- 
ence of magnetic fields. At distant cosmological objects the redshift may move the energy 
down into the optical or even the synchrotron frequency range and absorption of radiation 
on the line of sight by interfering material may suppress part of the spectrum. 

2.1 Pulsar Wind Termination Shocks 

Pulsar winds have probably the highest (or at least with certainty very high) Lorentz factors 
of the order of up to F ~ 10*. Compared with the energy of this flow, the wind itself is 
in good approximation cold. In addition pulsar winds are strongly magnetised with ratio of 
magnetic to total kinetic energy a ~ 10^^-10". Pulsar winds are expected to consist of 
electron-positron pairs (e+,e^), although ions may also be mixed into the winds, in case 
they are coUimated, particularly at their edges where ions may enter the wind from the the 
external medium by some mixing processes. 

In most cases (like the Crab pulsar, for instance; see [Gallant et aipOOS) ! the rapidly 
rotating pulsar of angular frequency Q, rotation period P = In/Q, moment of inertia Ip, 
and spin-down luminosity Lp = IpQQ which causes its windy outflow is embedded into 
the comparably slow expansion of the supernova which created the pulsar and into which 
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the pulsar is immersed. Hence, the pulsar wind environment is the supernova wind. The 
interaction between the pulsar wind and the upstream Supernova Remnant (SNR) material 
causes a pulsar-wind bow-shock system, which is located outside the outer boundary of the 
wind. Behind the bow shock one finds some turbulent (downstream) transition region which 
is probably terminated by a contact discontinuity (or, if magnetised, a tangential discontinu- 
ity). Inside of the latter the pulsar wind develops a termination shock that, in the frame of 
the wind, moves inward toward the central pulsar. Both these shocks, the bow shock and the 
termination shock, are believed to be relativistic and collisionless. 

The expected physics of the pulsar wind has been described by [Kennel and Coroniti| 
(|1984a') who, in analogy to the non-relativistic shocks in the heliosphere, assumed that the 
cold pulsar wind thermalises at the termination shock (see the sketch in Figure [T}. The 
plasma downstream of the termination shock is hot. Its sound velocity Cj/c w \/3 is a frac- 
tion of the speed of light, and its luminosity at the radial location rjs of the termination 
shock is Lts — 4;rF'A^ir^^m,c^(l -|- CTi). It is believed that the magnetisation ratio CTi = 
B^/fiofiiNic^r ^ 1 is small, where A'l is the upstream density, B magnetic field, and m; the 
ion (proton) mass. This yields a downstream pressure 3^ts ~ \r^N\miC^ = Ljs/^Ttcrjg. 
Simple arguments then suggest ( [van der Swaluw et al|2003| l that this pressure is constant up 
to the bow shock. In principle, however, as long as the pulsar has not broken out of the SNR, 
the evolution of the entire system depends in the external SNR-medium which determines 
the conditions at the bow shock. Recently, the spectral evolution of this internal medium has 
been investigated byl Bucciantini et al| ( |2011[ l who found an injection spectrum following a 
broken power law reproducing very well the observed radiation. After break-out from the 
SNR into the interstellar medium the dynamics and shock properties are determined by the 
conditions in the interstellar gas. 



2.2 Very-High Energy Gamma Rays 

Another indication for the possible presence of relativistic shocks is the observation of 
Gamma Rays in the very-high energy range <£' > TeV. Radiation of this kind has been ob- 
served in different manifestations in the Galaxy ( [Aharonian and et al.|2006c| ). The origins 
of diffuse TeV emission in its central region are probably dense molecular clouds which 
stand in the way of high-energy Cosmic Rays that most probably have been accelerated by 
shocks. The most probable three galactic sources of high-energy Cosmic Rays are SNRs, 
Pulsar Wind Nebulae (PWNe) and binary systems (BSs) (for a recent of detection > 100 
GeV gamma rays from binary pulsars cf. [Aharonian and et al.|2009bl ) which may produce 
the required ultra-high energy Cosmic Rays. Extragalactic sources can be active galactic nu- 
clei (AGNs and Blazars) ( jAcciari and et al.|2009| l, radio and starburst galaxies (|Acero and 



et al. 12009) 1 and, of course, clusters of galaxies (examples are given in [Aharonian and et al. 
2009a|c I. All arguments rely on the diffusive shock acceleration mechanism (DSA) in its 



various more or less sophisticated versions (cf., e.g. Bykov and Toptygin,1985)[Blandford| 



[and Eichler|'T987l[Malkov and Dmry|200T][Dmry et al|2001[ and others). 



Galactic sources 

Approximately 10% of the total mechanical energy release in the Galaxy stems from SNRs. 
This suffices to generate the high-energy Cosmic Rays at the obser ved energy density o f 



10 ^ J m . For an example observed recently with Chandra see Bamba et al (2003 



Young SNRs, in particular, are strong gamma-ray emitters both by the hadronic Cosmic 
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Fig. 2 Left: gamma-ray image above 300 GeV obtained witli the HESS telescope array in 2004. this image 
shows the young SNR source RX J17. 13.7-3946 which was observed in China as a 'guest star' in 393 AD. 
Colour coding is linear based on photon counts. In the lower left lower corner of this picture, which was 
adapted from Aharonian and et al. 1 2006a I the point spread function is indicated. The high energy gamma 
emission shows an asymmetric shell structure that is typical for the interaction of the SNR blast shock wave 
with an external molecular cloud. Right: Interestingly the spectrum is very hard, and is not a simple power 
law; moreover, it shows no angular dependence. A tentative model explanation has been given by assuming 
a modification of the diffuse shock acceleration theory jMalkov et al|2005) . 



Ray component as by 10 TeV electrons, which emit synchrotron or inverse Compton 
radiation and are thus strong candidates for relativistic shock acceleration. An example of 
the gamma radiation from a young (i.e. a historical) SNR is shown in Figure[2] 

Terminating pulsar winds are the other connection between galactic very-high energy 
gamma rays and shocks. Pulsar winds are ultra-relativistic. Ending up in the diffuse Pulsar 
Wind Nebula they generate a termination shock whereby electrons will probably be accel- 
erated up to 100 TeV energy or even beyond. These electrons should produce TeV gamma 
rays when undergoing inverse Compton scattering on the diffuse radiation fields and push 
the radiation energy up into the very-high energy gamma range. Clearly, inverse Compton 
scattering by the abundant Cosmic Microwave Background photons will also create lower 
energy radiation below 100 GeV. Radiation of this kind has been observed from the Crab. 
Its spectrum is unique in the sense that it extends over 21 decades in frequency or energy, 
from radio frequencies to very-high energy gamma rays implying that electrons are accel- 
erated up to > 10"' TeV. In a stochastic diffuse acceleration process they assume random 
pitch angles so that the bulk of the electron energy is released as radiation in the strong 
pulsar nebula magnetic field of ~ lOOjUG as synchrotron emission. At the same time inverse 
Compton scattering generates the very-high energy gamma emission. Besides the Crab neb- 
ula, which is exceptional in its radiation, CANGAROO and HESS have identified four other 
Pulsar Wind Nebulae which are also efficient gamma emitters: Vela X ( Aharonian and et al. 
2006e| l, MSH 15-51 ( |Aharonian and et al.|2005ai ), the PSR J1826-1334 nebula (Aharonian 
and et al.|2006d| l, and Kookaburra ( [Aharonian and et al.|2006bl l. The radiation source in all 
these cases is clearly separated from the corresponding pulsars indicating that the emission 
does not come from the pulsar itself but rather from the nebula caused by the deceleration 
of the pulsar wind at the remote location of the reverse termination shock as was suggested 
in the simple pulsar wind MHD model of |Kennel and Coroniti| ( [T984a|b[ l. 
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In pulsar nebulae the particles are accelerated in the external shock region, but in com- 
pact binaries one expects complete thermalisation of the plasma. Yet, the observation of 
very-high energy gamma rays from such objects as well suggests that this assumption might 
not be true. As long as one sticks to the shock acceleration picture the resolution of this 
puzzle can be seen in a number of internal shocks located inside the jet and close to the 
compact object which serve as efficient accelerators for particles. This can be the case when 
the object does, indeed, eject jets. Such objects, called micro-quasars, have been found in 
the galaxy, too (Mirabel and Rodriguez [T994| ). The binary object LS 1+61 303 ( [Albert and| 
|et al.|2006) ) belongs to this category of objects, it has radio emitting jets located at a distance 
of w2 kpc. It emits gamma rays up to energies of 4 TeV, possibly indicating that external 
termination (or internal) shocks are created. The gamma-ray luminosity of this source above 
w200 GeV is sufficiently high at w 10^^ J s^' to support this assumption. 

Finally, of the galactic sources the central galactic engine near Sgr A*, a probably super- 
massive black hole and its environment, is of interest. TeV gamma radiation with no time 
variability has been observed from this region where, however, a pulsar is also in closer 
vicinity, where ejection of material could cause shock generation in jets and termination 
shocks ( [Aharonian and Neronov|2005b|a| l. 

Extragalactic sources 

As mentioned above, potential extragalactic candidates for relativistic shock signatures in 
very-high energy gamma rays are AGNs, radio and starburst galaxies, and galaxy clusters. 

Of the AGN candidates, AGNs with relativistic jets close to the line of sight (with F 
few times 10 and a ~ 10^^ or less), so-called Blazars and among them the BL Lac objects, 
are the strongest emitters of observed radiation^This is believed to be caused by relativistic 
Doppler-boosting due to the line-of-sight orientation of the jets. Emission in these objects 
covers the entire range from radio to gamma rays. TeV gamma rays have been detected from 
the BL Lac objects. Their radiation is highly variable in brightness and polarisation which is 
a signature of the presence of ultrarelativistic electrons, which must have been accelerated to 
these energies. Variability in gamma flares in some of these objects has been found to be as 
fast as several minutes (for instance in Mkn 501, as reported by |Albert and et al.|2007[ l. In the 
flaring state the gamma-ray spectra of these objects are generally very hard. They can well 
be described by an exponentially truncated power law ANy/dS cx S^" e.Tup[—S J Sq) with 
power a ~ 2 and truncation energy Sq. This is shown in Figure [3] Such spectra are typical 
for shock-accelerated particles causing the radiation. More distant BL Lacs seem to have 
steeper photon spectra probably caused by intergalactic absorption. However, generally, the 
steepness seems to be intrinsic, allowing for hard spectra but not harder than o; = | . From 
energy requirements and statistical mechanical theory such a power-law index is an absolute 
lower for any phase space distribution. 

One principal point to be noticed in all these observations is that interpretation of the 
gamma emission by shock-accelerated particles encounters serious difficulties when we as- 
sume that the spectra are caused by hadrons, i.e. from protons or other baryons. The emis- 
sion must be caused mostly by leptons (electrons and positrons), and this poses the problem 
of how these particles are efficiently accelerated to the observed extremely high energies 

' We may note at this point that it was the observational evidence from relativistic explosions in AGNs that 
led |Blandford and McKee|jl976^ to develop their ultra-relativistic (coUisional) self-similar (non-magnetised) 
fluid model of blast-wave driven ultra-relativistic shocks on which most of the subsequent literature was 
based. 
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Fig. 3 Left: gamma-ray spectra for Mkn 501 and Mkn 421 during gamma flare states. Such spectra can be 
composed of exponentially truncated power laws (adapted from |Aharonian and et al.|2002^ . The lower part 
of the figure shows the ratio of the emission spectra of the two objects. Right: The spectral cut-off energy as 
function of gamma photon irradiation for Mkn 421 showing that the irradiance F{S >2) TeV increases with 
truncation energy S^ut = S'o according to a weak power law (data taken from |Aharonian and et al.|2005b^ . 

by shocks. This is a general problem which has not yet been solved satisfactorily. It re- 
quires either direct leptonic acceleration by shocks, which is a problem in itself, and pos- 
sibly might imply a coupling with other acceleration mechanisms until shock acceleration 
may set on for leptons, or it requires an efficient redistribution of energy from high energy 
shock-accelerated hadrons to leptons. 

Other candidates are radiogalaxies. TeV gamma rays have been detected from the nu- 
cleus of the nearby galaxy Centaums A and from M87 (which is 16 Mpc away in the Virgo 
cluster). In M87 the gamma rays may be produced either in the assumed central supermas- 
sive black hole or a hot knot (HST-1) roughly 100 pc displaced from the nucleus. The latter 
possibility is not improbable as the variability of the gamma emission correlates with that 
of the hot knot. Since hot knots are believed to be related to the termination of jets flowing 
out of the central engine, relativistic shocks would naturally be involved in the generation of 
this radiation, again implying violent particle acceleration. 

In starburst galaxies like M82 and NGC 253 there is plenty reason for strong winds, 
shock formation, particle acceleration, and generation of very-high energy gamma emission 
simply because of the high expected supernova rate. Indeed, gamma rays up to TeV energy 
have been detected from these objects. Further candidate objects are ultra-luminous infrared 
galaxies like Arp220 (located at w 70 Mpc). Finally, galaxy clusters containing large num- 
bers of either of these objects are clearly candidates of relativistic shocks. Relativistic shocks 
in those clusters are, however, probably not to be found in the intracluster medium, they are 
expected rather to belong to the galaxies which are members of the cluster. Observation of 
such shocks is however hindered by the comparably high thermal X-ray flux from clusters 
and the high absorptivity of the intracluster plasma which might prevent the observation of 
gamma rays. 
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Fig. 4 A so-called 'Fireball' sce- 
nario in which a GRB is caused 
by the relativistic outflow from the 
merger of two massive objects (af- 
ter [Waxman 2006). The relativis- 
tic Lorentz factor is assumed to 
be of order F ~ 300. Gamma rays 
may emerge with high luminosi- 
ties of the order of Ly ~ lO''^ J 
s"' from internal shocks in the 
highly relativistic outflow as syn- 
chrotron emission from electrons 
at >MeV photon energy. When 
the flow interacts with the envi- 
ronment, an external collisionless 
shock is generated which causes 
the synchrotron GRB afterglow. 



2.3 Gamma Ray Bursts 



Gamma Ray Bursts (GRBs) are a particular class of gamma-ray sources. We consider them 
here separately from the above-mentioned gamma-ray evidences for relativistic shocks of 
Lorentz factors F > 10^ with grossly unknown magnetisation ratios a. As their name tells, 
they are brief gamma-ray splashes that last only for short times (minutes to days) and are sin- 
gular events, probably of extra-galactic origin. Since their discovery in the sixties and their 
identification with remote events in the eighties, GRBs have been in the focus of astrophysi- 
cal interest partly because of their enormous brightness and luminosity, their time variability, 
nearly homogeneous distribution over the sky, and partly because of the suspicion that they 
would have cosmological relevance. Meanwhile a large amount of review articles on GRBs 



has accumulated in the scientific literature (as for a selection cf., e.g.. 


Gruzinov and Wax-| 


|man| 1999a| |Piran| 1999b| |2004| |2005l |Meszaros|2002| |2006l |Waxman 


20061 |Nakar|2010|(. 



Gamma Ray Bursts are believed to be emitted when stellar-mass black holes form in the 
Galaxy (or in other galaxies as well) by accreting stellar mass objects. Birth of such black 
holes is accompanied by relativistic outflows for which numerous observational evidence 
has accumulated from interstellar radio scintillations ( |Frail et al|1997[ l, apparently superlu- 
minous motions in afterglows ( [Taylor et al|2004| l, and afterglow light curves ( jHarrison et al| 
[T999|. 

Though it still remains unclear what the real source of GRBs ultimately is, the most 
probable models of generation of Gamma Ray Bursts are based on the assumption of ultra- 
relativistic collisionless shocks ( |Waxman]|2006| l that should be produced at least at some 
later stage in highly collimated ultra-relativistic outflows. From that point of view. Gamma 
Ray Bursts offer themselves as the ideal test case for ultra-relativistic shocks, their effects 
and properties. This view has been strongly pushed in the literature. Here we briefly re- 
view the supporting arguments. We will later, in the discussion of the theory of collisionless 
shocks, return to GRBs as a field of application of ultra-relativistic collisionless shock the- 
ory which is mainly based on numerical particle-in-cell (PIC) simulation techniques (for an 
early comprehensive review cf. |Dawson|| 1983) 1 which today can be applied to quite large 
numbers of particles per cell, extended boxes, up to three spatial dimensions and even inclu- 
sion of realistic ion-to-electron mass ratios = nii/mg = 1840. 
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Progenitor, central engine, outflows, jets, radiation 

Observationally GRBs can be divided into long (L) and short (S) GRBs. Only the former 
have been shown to be related to supernovae, exploding within ~1 day of the supernova 
(SN). Such a connection was ruled out for the latter. It is thus suggested that LGRBs occur 
in host galaxies with high star formation rate. Establishing such a connection immediately 
implies that LGRBs are at least accompanied by external shocks which, like in SNs, are 
generated when the outflow blast interacts with the interstellar medium in the host galaxy. 
In addition, internal shocks contained in the outflow may also exist, but a large part of 
particle acceleration and radiation in the observed spectrum is probably generated in the ex- 
ternal shock. Estimates of the relativistic Lorentz factor of the outflow univocally agree that 
F > 30 if not r ^ 100 as suggested by the detection of GeV gamma-ray photons. These 
values certainly correspond to angularly highly collimated (i.e. beamed) outflow as other- 
wise the released energy would exceed 10^^ J. Beaming has been confirmed by afterglow 
observations of GRB 030329. 

What concerns the location of the prompt emission of radiation, it seems to be certain 
that it comes from a radial interval of (10^ cm to lO'^m) < r < (lO'^cmto lO'^m) away 
from the central object, converting > 10% of the kinetic energy of the outflow into radia- 
tion, most of it into sub-MeV gamma radiation. This range suggests that it is emitted from 
the outflow itself, probably from internal processes as, for example, particles accelerated by 
internal shocks and emitting synchrotron radiation. This seems to be the case, because ob- 
servation has ruled out inverse Compton radiation, as this would require too strong magnetic 
fields inside the jet. 

Since the connection to SNs is missing for SRBs, their origin is tentatively attributed to 
the merger of either two neutron stars or a neutron star and a black hole. Also magnetars or 
proto-magnetar systems are in the focus. Under such circumstances high energy relativistic 
jets will be emitted from the central engine which would explain the high energy release 
of up to w 10^^ J (for the recent observation of the most distant cosmological GRB cf. 
[Chandra a nd et al.|20I0| l. All these assumptions are highly hypothetical so far; they could 
possibly be confirmed by presence or absence of gravitational wave signals emitted by them. 
However, in this case the generation of relativistic plasma shocks is less certain, even though 
the merger of the heavy objects is a violent process in the course of which blast waves might 
be released that propagate radially out, accelerate particles and terminate in the interstellar 
medium. 

Gamma Ray Burst Afterglows 

GRBs are followed by long lasting afterglows covering the entire radiation spectrum over 
eight orders of magnitude in energy or frequency, from X-ray energies to radio wavelengths. 
This radiation is almost certainly emitted in the interaction with the interstellar or galac- 
tic medium of density A' which terminates the outflow blast from the central object (for a 
more recent account cf., e.g., |Waxman|2006{|Katz et al|2007[ l. This interaction results in the 
generation of a pair of shocks: A strong external forward shock of radius r^ and a mildly rel- 
ativistic reverse shock (the equivalent to a termination shock in stellar wind outflows) that, 
in the frame of the forward shock, propagates inward. The reverse shock is, however, short 
lived in this highly bursty or time-dependent case. 

One believes that within minutes of the outburst the shock Lorentz factor becomes of 
the order of F ~ several 10^ while decelerating to F ~ few 10 within one day, and to F ~ 
0(1) within few months, by then becoming non-relativistic. 
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It is believed, moreover, that upstream of the shock the magnetic energy density is neg- 
ligibly small compared with the kinetic energy of the flow, its upstream CT-factor (neglecting 
thermal energy density) probably being as small as (Ji = /IfioNiirijC^r w 10^'' to 10^^. 
If true, this would make the shock about non-magnetic upstream. Downstream the mag- 
netic field is that of the ambient plasma 5 ~ /xG and should be tangential to the shock. 
Since this causes problems with the condition of magnetic field continuity, such shocks nec- 
essarily must generate their own selfconsistent fields close upstream to the shock as well. 
Moreover, from the observation of the emitted synchrotron radiation, a is believed to in- 
crease across the shock in the transition to downstream by orders of magnitude to values 
G2 w \Q-^to G2 w 10^'. This increase is traced back to some - still badly understood - 
mechanism of generation of magnetic fields in the shock environment. 

For a while, usually a few hours after the GRB, these shocks accelerate particles until the 
particles start radiating. Thereby the external shock decelerates (assuming proper number 
~ AFN and energy ~ AF^NmpC^ density conservation of the freshly shocked plasma in the 
plasma rest frame). The shocked plasma is contained in a comparably thin shell of width 

The simple self-similar canonical ultra-relativistic blast-wave shock model ( [Blandford] 
|and McKeej 1 976 1) predicts a decrease of the Lorentz factor with shock radius r, in the shock 
frame given approximately by 

Distant observes should see the shock expanding as rs{t) ~ \Q^'^{tS'p/N)^ m where the 
proton energy Sp is measured in IC** J, the ambient density is in m^-', and the time t is mea- 
sured in seconds. This seems not to contradict observation of the afterglow radiation when 
assuming that electrons couple energetically to the protons (baryons) and obey a power law 
distribution dNe/de °^ with p > 2. Afterwards, radiation causes cooling the electrons 
while de-coupling them from the baryons. 

However, even though this might be true in general, the details of the afterglow do often 
not allow for the application of the simple blast-shock model. This suggests that not only the 
mechanism of GRBs is not sufficiently understood today, but also the relativistic or ultra- 
relativistic shock theory is not yet capable of describing the afterglow satisfactorily (see, 
e.g., |Nakar|2007||20T0| ). 

Recently it has been proposed (> Meszaros and Rees|2010[ l that Population III stars might 
be the most probable cause for GRBs involving external shocks. Population III stars are very 
massive first generation stars, which may collapse already at early times and large redshifts. 
Once they collapse they will produce relativistic jets and splashes of radiation which could 
be those which appear at GRBs. This proposal can be checked by inferring about the special 
properties ultra-relativistic shocks should exhibit when they are involved into the generation 
of GRB afterglow ( |Fox and Meszaros|2006]|Toma et al|20I0^ . 

It is almost certain that the afterglow radiation is emitted by electrons by the synchrotron 
mechanism. This again requires strong downstream magnetic fields which, by the observa- 
tions, should be strong enough to be close to energy equipartition on scales much larger than 
the nominal shock width which, for non-magnetic upstream flows is the upstream (relativis- 
tically invariant) ion inertial length A/ = c/cOpi w 100/-\/]V[ km, where A'l is measured in 
cm^-' and cOpi is the ion-plasma frequency. Observations then require that downstream mag- 
netic fields should be strong over downstream (index d) distances w IFct 2> ~ A,-, 
much larger than the shock width 4,, where r ~ 1 day, yielding ~ lO'^zi,, and electrons 
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should become accelerated by the shock to cosmic ray energies. This implies the clarifica- 
tion of downstream shock electron acceleration to the observed high energies and the high 
inferred (near equipartition) magnetic fields. These, too, pose problems that so far cannot be 
considered to have been solved to any satisfaction. 



2.4 Ultra-High Energy Cosmic Rays 

The origin of Cosmic Rays is still insufficiently understood even today. The urgency of 
understanding it better can be illustrated when one notices that the mean cosmic ray energy 
density in the Universe is of the same order of magnitude as that of the cosmic microwave 
background: ucr ~ ucmb- This might be a coincidence, however, it is difficult to understand 
what the origin of this agreement is. 

At least part of the Cosmic Rays is believed to be accelerated in relativistic or ultra- 
relativistic shocks in many objects, including the sources of GRBs and starburst galaxies. 
A recent HESS observation of Cosmic Rays from the latter (see |Acero and et al.|[2009| ) 
with very intense gamma-ray fluxes at 2.2 GeV, implies cosmic-ray densities three orders 
of magnitude larger than those emitted from the centre of our galaxy. At the same time 
such Cosmic Rays also contribute to the structure and dynamics of those shocks. Diffusive 



shock acceleration seems not to work at energies higher than 10^^ eV to 10'^ eV (Meli and 
|Biermann|2006[ ). Whether this limitation holds or not is unclear. 

Sub-relativistic Supernova shocks can accelerate Cosmic Rays presumably up to ener- 
gies of < lO'^ eV i Ginzburg and Syrovatskii|l969 Bykov and Toptygin 1985; Bland ford| 



|and Eichler|1987[ l. External shocks in jets of active galactic nuclei (Berezinsky 2008^ and 
LGRBs have been suggested early on (cf., e.g., |Waxman|1995a|b[|Vietri|1995[ Wick et al| 



|2004[ [Dermer and Atoyari||2006 ', and others) to be sources of (galactic as well as extra 
galactic) Cosmic Rays in the ultra-high energy range e ~ 10"* eV to 10^" eV. This might 
indeed be the case with the exception of the highest cosmic ray energies > 10^*^ eV for 
which no reasonable shock acceleration mechanism is known (cf., however, the - debatable 
- discussion in |Dermer|2006[ where it was argued that external LGRB shocks could as well 
generate particles at such energies). The limiting energy for shock accelerated protons has 
been estimated as 

'100 



where Lp is the undisturbed luminosity (in J) of the upstream flow and CTi = B\/2iiqw\ the 
ratio of upstream magnetic field energy density to the total energy density w\ of the flow 
including rest energy ( jWaxmanj 1 995a| l. Clearly, for the small CTi <S 0. 1 used in the previous 
section, one never reaches the wanted limit of lO^'' eV. 

Short GRBs can also accelerate protons to such high energies though at much reduced 
flux by internal-shock acceleration in the upstream flow, if the Lorentz factor s range between 



10^ lij 10"' and the magnetisation is at least moderate (Nakar 2007 i. For a general 
discussion of how to extract information on a retarding expanding external relativistic shock 
that can be extracted of a cosmic-ray spectrum one may consult a recent paper by jKatz et al| 
(2010 ). 

In the Galaxy there is another relation for the energy density of Cosmic Rays: ucr ~ 
/Ifjf) ~ ^PgVnirb' ^^^^ Pg density, and V,urb the average turbulent velocity of galactic 
matter. This relation says that the cosmic-ray energy density is comparable to both the mag- 
netic and turbulent gas energy densities - a strong indication that neither the magnetic field 
nor the gas can be considered independent of the presence of Cosmic Rays. 
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Conversely this mutual dependence also indicates that Cosmic Rays are in energy ex- 
change with plasma and magnetic field, and as well as shocks. In particular relativistic and 
ultra-relativistic shocks are the most probable cosmic -ray factories and are thus responsible 
for this energy exchange - in addition to general turbulence. 



3 Advances in Relativistic Shock Theory 

Almost until the start of the third millennium, reference to relativistic shocks has been based 
solely on the fluid approach as given by |BIandford and McKee] ( |T976l l, pc"Kee and Ostriker] 
( [1977| ) and others. The substantially more advanced non-relativistic shock theory and obser- 
vation in space physics l |Tsurutani and Stone|1985| l was mentioned only occasionally, albeit 
with the clear notion of not being applicable to astrophysical problems. 

The change came with the realisation that shock acceleration stagnated as long as ki- 
netic effects and magnetic fields were not included into theory. Including these effects was 
first proposed by |Gruzinov and Waxman| ( |l999b[ l,|Medvedev and Loeb (1999) and' Brainerd| 
pOOOj ) to account for the generation of the large magnetic fields which are required in the 
shock-radiation models of Gamma Ray Bursts ( [Piran|I999b|a| ). In the following we briefly 
review the basic physical background of shock theory and the current state of the art with 
the advances in the latter coming mainly from numerical simulation. 



3.1 Fluid Dynamics of Relativistic Shocks 

As noted above the study of relativistic and ultra-relativistic shocks began with the investi- 
gation of non-magnetised fluid dynamic shocks ( ,Blandford and McKee|1976l|1977[|McKee| 
|and Ostriker|1977[|Blandford and Ostriker|1980| ). Since it became clear that there are no non- 
magnetic shocks in the Universe, at least no visible or otherwise detectable non-magnetic 
shocks, we do not review these seminal works but go straight to a brief consideration of 
relativistic magneto-fluid dynamic (magnetohydrodynamic) shocks. 

Relativistic Ideal MHD Rankine-Hugoniot Conditions 

Shock jump (Rankine-Hugoniot) conditions in one-dimensional ultra-relativistic £ ^ {m + 
pY)c^ non-magnetic (hydrodynamic) fluid flows is the volume, S the energy) for ap- 
plication to blast shock waves that may presumably be generated in radial outflows from 
astrophysical explosions were given by Bland ford and McKee| ( |197^ in their most simple 
form and have been widely applied. Their extension to spherical relativistic ideal MHD flows 
is due to |Emmering and Chevalier] l |1987| l. 

In terms of the velocity four- vector = (u' ,1/) = (l,jS''),;S = V/c, with rest-mass 
density p = Nm, proper number density A^, average rest mass per particle m, magnetic field 
H = Bq/^/Atz, with magnetic field So and gas pressure 3^ in the fluid frame, enthalpy w = 
pc^+g3^, where as before g = f/(f — 1 ) is the ratio of the adiabatic index y, the temperature 
T = £y'/p, and putting c = 1. Working in the shock frame (index s) where the quantities are 
assumed stationary, the jump conditions (written in one line only) become 

[[piQ = [[{w + H^)u'X]] = [[(w + H^K^ + ^+ = [[^/p]] = (3) 
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The double-brackets [[]] indicate the mismatches of quantities to both sides of the shock: 
[[A]] = A2 — Ai where the indices 1 and 2 refer to upstream and downstream of the shock. In 
the observer's (fixed shock) frame the electromagnetic fields are given by B = m'Bq and E = 
i/Bq. Any normal (in this case radial) magnetic field component remains continuous because 
of the vanishing divergence condition. In the non-magnetic limit H ^0 the jump conditions 
become the hydrodynamic jump condition^ and in the ultra-relativistic radially expanding 
symmetric case of very large bulk Lorentz factor F S> 1 are just another generalised implicit 
representation of those given by ,Blandford and McKee| l (T976p . 



Magnetised Plane Perpendicular Shock Solution 



For a magnetised, plane and strictly perpendicular shock a set of explicit shock relations can 
be derived from Eq. ^ best when working in the downstream frame of reference (cf., e.g., 
|Amato and Arons^20 06). The magnetic field B = Sz points in direction z along the shock 
(A:,z)-plane, the shock moves at velocity <C jS^/, = Vsh/c < 1 upstream in direction —x, 
and the upstream flow is cold and has Lorentz factor F. As before with magnetisation ratio 
a = B\/ jM.)imc^N\r , m = Lj/Mj, and hot downstream plasma of ultra-relativistic enthalpy 
W2 w ^27/ (7— 1)> one finds 



N2 B2 1 + A 



'sh 



(= Rs, the shock compression ratio) (4) 



mc'r{l + G)Ni = 3^2 + ^ = j^r(^ + ^] (5) 



2^0 V7-1 2^0 

Here, the ratio of specific heats is 7 = | in ultra-relativistic isotropic and 7=5 two- 
dimensional cases, respectively. (These equations are not applicable in the strictly ultra- 
relativistic case for a shock with jS.,/, :^ 1 , as this would require 7 < 5 •) After elimination of 
P2 and B2 a quadratic equation is obtained for the shock velocity jSj/, and its solution can be 
used to find the normalised downstream temperature T2 



(6) 



iS^/' = ^Y^|(f(l + 5^)-l) + [(7(l + ^cT)-l)'+4a(l+cT)(l-if) 

^^rA.(i-fl^) (7) 

(with temperature in energy units) showing that in this simple case the step across the shock 
is determined for all quantities by three parameters only: a, 7 and F. Since T2 > 0, this also 
implies that 1 > jS,/, > o/{2 + a), a condition that is easily satisfied for weakly magnetised 
upstream flows. 



- One should cautiously note that physically the non-magnetic (hydrodynamic) and magnetic (magneto- 
hydrodynamic) cases are fundamentally different. In the latter case, the gyro-magnetic motion of the charged 
fluid particles causes particle correlations even in the absence of any binary collisions thus justifying the 
coUisionless ideal MHD assumption. Such correlations are missing in the hydrodynamic case which chal- 
lenges the coUisionless assumption, i.e. implies the presence of some kind of collisions. Moreover, applying 
the non-magnetic approximation to plasmas consisting of charged particles is particularly interesting. Un- 
der certain conditions (see below), such hypothetic non-magnetic plasmas act self-magnetising on very short 
time scales, a posteriori justifying the magnetised-shock approach. However, the validity of MHD may be 
questioned in principle because the large population of shock-accelerated particles around relativistic shocks 
(cf., e.g., |Gruzinov|2001| requires a kinetic treatment. 



Collisionless Relativistic Shocks 



15 



The Case of Spherical Symmetry 

Though a sufficiently thin shock transition can, locally, always be taken as plane, in the cases 
of rapidly expanding pulsar nebulae (or also ultra-relativistic jets) spherical symmetry has 
traditionally been assumed to be more appropriate ( Kennel and Coroniti,1984a|b,,Emmering| 
[and Chevalier|198"7] l. In spherical coordinates, the last of conditions ([3]l holds automatically 
because the radial mass continuity and azimuthal magnetic field equations assume identical 
form 

di{r-^Hu') + r-^dr{r^{r-^Hy} = 

suggesting that r^^H varies like p and thus the ratio H/rp = const is conserved everywhere 
(also holding across the shock). These spherical jump conditions have been widely discussed 
analytically(cf., e.g.,I Kennel and Coroniti|1984a|b[ who considered a stationary shock) and 
have also been solved numerically for radially diverging flow under various assumptions 
imposed on the shock-upstream parameters (cf., e.g., |Emmering and Chevalier|1987| l. 

Traditionally two different cases have been considered corresponding to Vsh > V\ and 
Vi/i < Vi , where Vsh is the shock speed, and Vi is the shock-upstream velocity. The former 
case describes so-called 'forward', the latter 'reverse' shocks ( [Blandford and McKee|1976| 
[Kennel and Coromti|1984a]bl l. In the first case the shocked fluid is behind the shock, in the 
second case it is ahead of the shock, while in both cases the shock moves outward towards 
increasing radii r. 



3.2 Self-similar shock structure 



Self-similarity at ultra-relativi stic non-ma gnetic astroph ysic al blast waves w as introduced by 

Below we first 



Blandford and McKee 


1 1976 i, following 


Sedov 


\959\ and 


Taylor ( 


1950P 



refer to an approximate self-similar treatment of the radially symmetric shock structure, and 
then treat the general self-similar case. 



Approximate Self-similar Solutions 

The radially symmetric, self-similar equations contain a singularity that is located at the 
critical velocity ratio rjc = Vc/Vsh = jic/jish with Vc the critical velocity. At the shock T) = 1, 
and hence the size of the shocked region is limited between the location of the shock and the 
location of the singularity. For reverse shocks one has T)c — 1 < {Pshr^^i^ + As/i)} This 
yields for an ultra-relativistic shock with jSj^ ~ 1 and 2> 1 that the reverse shock-shocked 
region is very narrow: 

T],-i<ir,-2 (8) 



Because of this narrowness the value of Vc is practically identical to the velocity 152s in the 
shocked region adjacent to the shock. Thus in the fixed frame (after Lorentz transformation) 
one has for the width of the shocked region 



\ric-l\-p2s/Vshrs^{l±VshP2s) 



(9) 



It may be worth mentioning that Sir Geoffrey Taylor's paper on self-similarity solutions of blast waves 
was originally written as a classified paper during the Second World War in 1941. 
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Fig. 5 The postshock pressure P2 (thick lines) and magnetic field B2 (thin lines) determined from the shock 
jump conditions (after |Emmering and Chevalier|1987l for a typical SNR like the Crab for different values 
of a and in dependence of the radius measured in SNR radii Rq. The assumptions made are that ijj ~ 11 
is in the extreme relativistic range, and the velocity at rj^, the boundary of the shocked region, is 2000 



kms ' corresponding the Q = 6.7 x 
f = tanh"'(H'7M'),?7 = r/fj = r/Vsht. 



10 . The position of the shock front is taken at 21}^ pc. Here 



The minus sign holds for forward shocks in this expression. The value of jS2j follows from 
solving the jump conditions. For small a (weak magnetic fields) the shocked-region width 
depends very weakly on a. The interesting case however is that of strong magnetic fields 
(large a) where jS2j ~ 1 — 1/2ct. This yields for forward and reverse shocks, respectively, 
the scalings 

l-T?c./or = 2a((J+r/)-\ T],,,,„-l~ ir,-2 (10) 

Hence, for reverse shocks the width shrinks with increasing F,. For forward shocks, 1 — rjcj'or 
grows with a and becomes 1 for ct ~ F^. Since forF, 1 one has /32.5 ~ (rj^ — o)/{r^ + a), 
the immediate post-shock velocity becomes small (non-relativistic), so that T]c.for = 0, which 
implies that similar to the non-relativistic case the entire region behind the shock contains 
shocked plasma. This is in contrast to reverse shocks where the width of the shocked plasma 
region shrinks as r]c.rev 1- Moreover, no forward shocks exist for such strong magnetic 
fields that a > F^, - at least in the relativistic MHD picture. 

The behaviour of shocked density N2, pressure P2, and shock-tangential (azimuthal) 
magnetic field B2 are of particular interest. For weak magnetic fields, i.e. small cr, the density 
in the shocked region is about constant. An application of the self-similar jump conditions to 
the expanding shocked matter in the Crab nebula has been given by jEmmering and Chevalier] 
( [T987| . 

Figure[5]provides a synopsis of the evolution of the shocked magnetic field and pressure 
following from a simple Crab shock expansion model as function of distance and for a 
couple of values a. The pressure (thick lines) drops to zero at the edge of the nebula while 
the magnetic field (thin lines) increases with radius toward the edge. The model which fits 
the observations best is that with ct' ~ 16 which corresponds to a magnetic field of i? ~ 
2.6 X 10"'* G at the edge. 
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Self-similar Shock Structure: General Approach 

Recently, in view of an application to the afterglow of Gamma Ray Bursts, [Katz et al| ( |2007| l 
developed a fully self-similar theory of collisionless shocks including (self-generated) mag- 
netic fields. The basic assumption of this approach starts from the observation of afterglow 
radiation which, in order to be emitted as leptonic synchrotron radiation, requires a strong 
downstream magnetic field with energy density close to equipartition in a region of extension 
L^/ 2> 4j ~ A/. The large observationally suggested magnitude of = Xcorr, interpreting it 
as a shock-magnetic correlation length, is taken as evidence for self-similarity of the shock 
structure. This approach is based on a separation of the plasma particle distribution into 
two self-similar components, one non-thermal, describing the shock-accelerated high en- 
ergy particle component (reflected particles, accelerated particles, a cloud of Cosmic Rays 
etc.), the other group are 'thermal' particles belonging to the flow component and are treated 
as a (magneto-hydrodynamic, in the end ideal) fluid component of vanishing resistance. 

The dynamics of nonthermal particles of species s = {e^,i) is described by the particle 
distribution functions Fs{t,p) which are governed by the set of (relativistic) Vlasov-Maxwell 
equations with current density J and single particle velocity v(p) respectively defined by 

J^l^eJ dV(p)^.(p), v(p) = '"^^ ^ (11) 

s J y'm^c + 

and far upstream (z < 0) boundary conditions 

(E,B)|,^_„ = 0, F,\,^^„=rN,iS^{p^)5{p,-rVun,) (12) 

A further and crucial assumption for the development of the self-similar theory is that 
the self-similar shock structure is strictly stationary, i.e. the non-stationary phase has settled 
into a stationary final state. This assumes that the time to reach a stationary state is much 
shorter than the time until the shock resolves into the environment by interaction with the 
ambient surrounding medium, in afterglows the time of roughly one month. The further 
evolution of the shock within this time window is considered to be self-similar and slow. 

Under these conditions, it is argued that, for small Gi = [Xi/rdjh)^ < 10^^, the upstream 
magnetic field plays no role anymore in the shock structure and the final well-developed self- 
similarity. This seems to be suggested by numerical simulations in pair (e^) plasmas (cf., 
e.g. , [Spitkovsky [2005 [ who reports simulations with Lorentz factor F ~ 30). The upstream 
field is just important in the initial shock formation process only|^ Under this assumption 
on the 'stationary shock dynamics', there will be no average large scale stationary magnetic 

field other than the root-mean-square magnetic field B{z) = E,i?,7(zixx = 0,At = 0,z)]2 
at normal distance z from the shock. This field results from the shock generated magnetic 
fluctuations and their non-zero field correlations 5,y. The average magnetic energy density 
at distance z is {S'siz)) = B^{z)/2^1q. 

** This claim is crucial as it enables to develop a stationary self-similar picture of the shock. It could, 
however, be questioned as it is by no means certain that a shock can reach a stationary state. The most sophis- 
ticated numerical PIC simulations of non-relativistic shocks, which do suffer much less from the problems 
encountered in simulating relativistic shocks, show that shocks are non-stationary on fairly long time scales 
of reformation (cf. Balog lTand Treumann|201 though it has also been suggested jLembege et al|2Q09) from 
non-relativistic two-dimensional simulations that reformation might be an artefact when restricting to one 
dimension only. Shocks also generate their own upstream magnetic field in various ways by instabilities in 
which the shock accelerated particle component is involved. This may imply that the initial upstream field is 
less important jKatz et ar[2007j , but in the stationary state a\ cannot be neglected anymore. In addition, one 
may expect that relativistic shocks do also experience quasi-periodic reformation cycles. If this is the case 
stationarity would apply to times much longer than the typical reformation time. 
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At large distances l^l A,- from the shock self-similarity implies that the correlation 
length Xcorr becomes the only relevant length in the system implying that the shock becomes 
scale invariant, if it starts diverging as Xcorr/h °°- This is assumed to be the case in 
self-similarity, yielding that all physical quantities become scale invariant power laws with 
particular power law indices aj according to ( |Katz et al|2007 1 

B(x,r) = ^) , F.(x,p,0 = ^«/, (^|, ^, ^) (13) 

where £, = ?icon-/^, and Ao is any fixed reference scale as, e.g., A,-. Denoting the non-thermal 
particle component by the index s = b, the problem then reduces to the determination of the 
mutual dependence of all the 'similarity indices'. For the thermal particle component one 
has for their gyro-radius the scaling rc = cp,i,/eB ~ ?i,corr and thus rdT^corr ~ "^corr ■ 
Consequently one finds that >«;,>— 1 for the bounds on the magnetic field scaling 
exponent. 

The remaining scaling exponents can be determined when inserting the above scaling 
relations for the field and distribution function into the Vlasov-Maxwell equation and requir- 
ing that the different terms should have similar scaling. Then one obtains for a, = 1, while 
the scaling power of the momentum should be ap = \ + aj,. This allows for the determina- 
tion of the relation between and 0!j from the scaling of the fast particle current density 
Jb as fixed from Ampere's law. The result is af + 3{ab + I) = oCb — I, or af = —2(2 + Ub). 
Since the average distribution function of the energetic particles in this case turns out to be 
a power law function, it can also be written in terms of the momentum as 

lim (F,(p,A,„„)) ^ "" = (P-) (15) 

p^"" \poJ \poJ 

where po is a reference momentum (see the discussion in |Katz et al||2007] l. Since ab < 
is restricted by the above arguments, the high energy particle momentum space distribution 
has power OCf/ap < —4, and the flattest particle distribution is obtained for = in a 
constant rms fluctuation field. Its correlation function {Bi{x)Bj{x + Ax)) and the spectrum 
of energetic particles that are accelerated by the self-similar coUisionless shock scale as 

{Bi{x)Bj(x + Ax)) ~ (Axf"", ^ ^ (16) 

asb 

with spatial distance x from the shock and particle energy £/,, respectively. Thus the self- 
similar shock-accelerated particle energy spectrum has power < —2. 

This self-similar theory is a very important step towards an understanding of relativistic 
coUisionless shocks. On the other hand, it is not clear whether the assumptions made can be 
justified other than by the heuristic arguments used by |Katz et al| ( |2007| l. These assumptions 
are the neglect of the effects of any initial upstream magnetic field on the final self-similar 
state in shock evolution and, more serious, the assumption that the self-consistently gen- 
erated magnetic fields can (i) indeed reach the required equi-partition amplitudes to which 
the self-similar theory refers, (ii) can fill the entire extended downstream and to some ex- 
tent also the region upstream of the shock. Clarification of these assumption is still far from 
being achievable. Nevertheless, the self-similar heuristic theory provides a tool for appli- 
cation to shocks. It makes a clear statement about the marginal power of the high energy 
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self-similar particle spectrum and it limits the self-similar exponent of the self-similar shock 
magnetic field at large distance from the shock. The range of this exponent is narrow. But it 
contains all the unknown physics the clarification of which is quite urgent and requires the 
investigation of the details of how a shock sustains itself and which processes are involved. 
Self-similar theory gives just the hint that the shock may sustain itself without providing any 
information about the underlying physical processes. 



3.3 Relativistic MHD Waves 

Shock waves form when the linear eigenmodes of the relativistic plasma evolve non-Iinearly. 
This happens whether Waves are excited by an external driver (a piston, the central explosion 
in the case of blast waves, etc.) or the shock results from wave steepening in high Mach 
number collisionless flows. In this respect, a relativistic plasma does not behave differently 
from its non-relativistic counter part. In order to infer about the nature of the shock and its 
evolution one thus needs to know the linear eigenmodes which are solutions of the linear 
dispersion relation. 

Observed astrophysical shock waves are large-scale. They evolve from large-scale low- 
frequency waves in plasma. One thus needs to find just the low-frequency wave branches of 
a relativistic plasma. The full theory requires a kinetic treatment based on the linearised rel- 
ativistic Vlasov equation. Observations suggest that magnetic fields B are involved, and one 
is mostly dealing with magnetised plasmas. To first approximation (and neglecting particle 
creation and annihilation) these may be treated by fluid theory based on the ideal relativistic 
MHD equation^ This is most elegantly done in the covariant formulation (cf., e.g 



Komis- 



|sarov|1999| . However, for transparency and applicational purposes the so-called 3-1-1 split 
formalism (cf., e.g., Keppens and Meliani 2008) is more appropriate because of its similarity 
to non-relativistic MHD. In this formulation one adds the ideal induction equation 

^-Vx(vxB)=0, V-B = (17) 
at 

to the set of relativistic hydrodynamic equations, which couples the magnetic field B to the 
dynamics via the velocity v in the Lorentz force E-l-vxB. One also adjusts the energy 
equation accordingly. 

Linearising around the uniform background quantities entropy ^o, average mass density 
po, and magnetic field Bq in the plasma rest frame (e.g., |Komissarov|1999j l one finds the rel- 
ativistic variants of Alfven and magnetosonic waves. In the 3h-1 split, assuming a polytropic 
equation of state with adiabatic index y, the linearised Fourier transformed equations for the 
first order quantities (index 1) with frequency O) and wave number k, read 

fflBi =Bok-vi-vik-Bo, k-Bi=0, Pi^ ^ k;V]_ ^^^^ 

Po (O 

= ^ (k+ (M + M+^(kBo.B,-B,k.Bo) (19) 

Here we defined the specific enthalpy h = c'- -\- gSoPg ' , w = po/i -|- Bq/hq, and g = 'y/(y~ 
1). The adiabatic index has the two limits y = 5/3 for non-relativistic and 7=4/3 for 

^ In the non-magnetic case one needs to refer to kinetic theory, however, in order to learn that magnetic 
fields can hardly be avoided, and to conclude that relativistic collisionless shocks should almost ever be 
accompanied by self-generated magnetic fields. This will be discussed in more detail below. 
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ultra-relativistic temperatures T, respectively (i.e. it depends on the internal microscopic 
velocities and not on the macroscopic speed of the flow). The solutions of these equations 
are compressible magnetosonic waves and the transverse Alfven wave vi ^ 0,Bi ^ 0,pi = 
= k ■ vi = (B, v)i - Bo = 0, with refraction index = k^c^ /aP' = jiQw/B^co&- 0, and thus 
Va/c = Bq cos 9 1 •^Hqw, where is the angle between wavenumber k and the undisturbed 
magnetic field Bq. 

Shock waves evolve from the compressible modes. Introducing the sound speed c\ = 
c^YPq ' ^0 1 h, their dispersion relation is 



CO - CO 



: (k-Bo) 
Mow 



Mow 



(20) 



We may note that for Wa c one has poh/w = 1 — V^/c^, and for propagation parallel to 
Bo, i.e. 0=0, the dispersion is the same as in the non-relativistic case. 

These expressions hold in the fluid frame. The fluid moves with relativistic speed V <c 
and bulk gamma-factor F > 1 with respect to the observer and shock frames. In these frames 
the above expressions for frequency, wavenumber and fields appear Lorentz transformed 
according to 



CO 



:r(co' + k'-V), k-k'=V 



coT r 

— ^ + 
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k'V 



B' 
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r ^r + 1 



r V B. 



-V (21) 



with r = I/i/l — V^/c^, and V is the relative velocity between the two (primed and un- 
primed) frames. For the stationary-frame Alfven velocity one obtains 



kc 



k v 



A.gr ■ 



k-V^ k-VA/r^fc 



(22) 



where y' A.gr is the Alfvenic group (wave energy flow) velocity in the stationary frame. 

Instead of transforming the magnetosonic wave dispersion relation into the observer's 
or shock frames it is more convenient to write it in terms of the phase speed \p = {co/k^)\i. 
which is a velocity and transforms (in terms of the transformation of the refraction index 
n = kc/ CO) according to 



1-2 _ 



j2 



r^{vp/c-]i-y/kcY 



r2(n-i-k-V/fc) 



c2 l-Kr2(Vp/c-k-V/fc)2-v2/c2 i+r2(«-i-k-v/fc) 



(23) 



In terms of the quantities in the moving frame this yields for the inverse refraction index of 
magnetosonic waves the quartic equation in the transformed frame 
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(24) 



Its solution can be given only numerically and depends heavily on the value of the Lorentz 
factor r of the fluid. An example if given in Figure [6] for the group velocity phase diagram 
of a fast mode in a medium at oblique relativistic propagation. There the group velocity 
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Fig. 6 The numerically calculated group ve- 
locity phase diagram (after [Keppens and| 
[Meliani 20081 for the fast mode in the ob- 
server frame, projected into the (vy, i'j)-plane. 
The cycloid line is the fast mode phase dia- 
gram; Units are normalised to po = 1 , Po = 1 
and Bo = 3. The velocity was taken as V = 
0.9c with direction angle 6 = 7t/4 (tip of ve- 
locity vector at blue dot). The wave phase is 
deformed by the relativistic motion. 



is shown in the projection(grey lines) of the three-dimensional wave-energy surface. Ob- 
viously, the transformation introduces the expected anisotropy of the wave pattern that the 
observer would see when having access to the magnetosonic phase and group velocities. 

As in the non-relativistic case, these fast waves may steepen and develop into a shock 
which is a thin interface between the relativistically fast super-magnetosonic unshocked 
upstream flow and the slow shocked downstream flow. Once this interface has reached its 
approximately stationary state the properties of the downstream flow can be inferred from 
the known properties of the upstream flow by using the appropriate jump conditions across 
the shock that must be satisfied. The following section gives an example of such conditions 
at a spherical shock with the flow being aligned radially with a radial magnetic field. 

3.4 Kinetic Wave Theory 

The fluid theory of the previous sections is not applicable to shock formation and shock 
structure. It just describes the jumps of the fluid parameters across the shock under ide- 
alised conditions: when the number of shock- reflected and shock- accelerated particles does 
not affect the shock structure. Already for high-Mach number nonrelativistic shocks this 
assumption becomes critical. For relativistic (and even more ultra-relativistic) shocks the 
efficiency of particle reflection and acceleration is very high, and the shock structure and 
jump conditions become modified by the presence of a dense energetic particle component 
which changes wave dispersion and affects the nonlinear evolution and steepening of plasma 
waves. Simple fluid theory ceases to be applicable then, and relativistic shocks would require 
a kinetic treatment. 

On the microscopic level the dynamics of collisionless plasmas follows from the Vlasov 
equation for the one-particle distribution function /(r,p,x) in the plasma frame 

a,/(f,p,x) + [jr,/(r,p,x]=0, |dV(r,P,x)=A^(f,x) (25) 

where is the single particle Hamiltonian, [a, b] is the Poisson bracket, and A' the particle 
number density to which the distribution function is normalised, p = myy is the relativistic 
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momentum of the particles, and 7(77) = ^1 + p^/m^c^ is the intrinsic Lorentz-factor of the 
plasma particles (which is a function of p and is not to be mixed with the external Lorentz- 
factor r of the bulk flow). 

The one-particle phase space distribut ion /(f,p,x) is Lorentz-invariant (see, e.g., Lan- 
dau and Lifshitz| 1975 [[van Kampen|l969p] as long as particle creation or annihilation are 



absent. Similarly, the pressure tensor P is an invariant as it is defined as the second differ- 
ential moment of /(p) in terms of the difference p — P between the single particle momenta 
and the bulk momentum P = mFV, which is the first moment of the moving distribution 
function being non-zero only in a streaming plasma. The (isotropic) scalar pressur^is given 
by P = ^ trace P (= NT for an ideal plasma). 

Since shocks are not in thermal equilibrium but require continuous driving by the up- 
stream flow and the interaction with kind of obstacles, they generate a large number of 
plasma waves in all frequency regimes from the plasma frequency (Ope = \J \Ke^'N j mg down 
to zero-frequency plasma modes. They also emit radiation at frequencies O) > COpe well 
above the plasma frequency by nonlinear processes and wave-particle interaction. Linear 
theory can be used to infer about the possible eigenmodes and their growth rates. However, 
the situation is complicated by the inhomogeneity and variability of the plasma in the shock 
making the application of the straight linear kinetic approach obsolete. The basic picture 
drawn by [Sagdeev] \ 1 966| l long ago for non-relativistic shocks is still valid in its basic lines 
and holds also for relativistic shocks what concerns shock formation and the responsibility 
of low frequency plasma modes. 

The complete spectrum of linear kinetic electromagnetic eigenmodes a)(k) = cOk excited 
in the shock and shock environment follows from the linear dispersion relation 



2 k^c^ 
it- = — =- = Det 



9 /I \ 

« p-+e(k,C0k) 



e(k,cOk) = l-F CT(k,cOk) (26) 



where I is the unit tensor, and e, ct are the dielectric and plasma conductivity tensors, respec- 
tively. The latter contain the entire plasma response and thus depend on the assumptions on 
the plasma model, including the various plasma components and sources of free energy. 

For purely transverse electromagnetic waves the longitudinal first term in the dispersion 
relation drops out, and the dispersion relation simplifies to just «^ = Dete(k,a)k). In this 
latter form it applies to the generation of the lowest frequency magnetic oscillations that 
constitute the shock, in particular to the description of the self-excited shock magnetic field 
by the |Weibel| ( [19 59')-zero-frequency magnetic instability that is claimed to take responsibil- 
ity for a much higher shock compression ratio Rs than provided by MHD shock theory. The 
latter applies to shock formation in the form of the |Fried| ( [1959] )-'filamentation' instability 
which describes the head-on interaction of two (collisionless) plasma beams, a model that 
applies to the spatially limited shock interaction region where the cold fast upstream flow 
mixes into the hot slow downstream plasma. 

The above expressions can also be written in covariant form (cf., e.g. |Dewar|1977[ l which 
in the relativistic case is sometimes more convenient for explicit calculation of the dielectric 
tensor expression. This makes use of the four-dimensional field tensor F^"^ = d^A^ — d^A^, 



^ The invariance of the distribution can be trivially seen if one recalls that the distribution is the probability 
of finding a particle in its phase space volume element. When correctly (covariantly) Lorentz transforming 
the latter from one frame into another frame this probability does of course not change. 

^ Thermally anisotropic but otherwise equilibrium plasmas require a different definition of temperature 
which is still debated upon. The coirect definition is analogous to Boltzmann's inverse definition and defines 
the inverse temperature as a four vector (cf. jvan Kampen|1968| . 
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with = (<^, A) the 4- vector potential, and the current 4-vectoi|^ In linear approxima- 
tion in Fourier-space k = = (a)/c,k) with scalar 4-vector product k-x = r]fivk^x^ = 
COf — k - x we have —iF^^{k) = k^A^ — k'^A^, and the connection between the current and 
potential of species s is given through An ^"^J^ = (■''A''v(^M^(^), where is the 

polarisation tensor contribution of species s which has the symmetry property k^ ^"^X^^ = 
^^''X^^kv = 0. The linear 4-current is the linear response of the plasma to a perturbation of 
the plasma distribution f{p,x,t), where p = = rrisUs is the 4-momentum, and u = u^= 
the 4-velocity, the overdot indicating differentiation with respect to proper time. This yields 
for the linear current '''/(A:) from the Vlasov equation 

^ ' J ' dp' ' \ ri^pk'^pP J 47: 

where co^ = ATte-Ns/nis is the square of the i-species plasma frequency, and fsi){p) the 
(known) undisturbed relativistic j-species distribution function. Thus the polarisation tensor 
becomes 

(■')A^W = -^/dV/,o(p) 

The set of equations for the linear 4-vector potential A' {k) then becomes D^' {k)Av{k) = 
with D^'ik) = {r\apk''kl^)r\f'' - k^'k" + X>''{k), where now = I, is the 

full polarisation tensor. 

The covariantly written but otherwise identical to the above non-covariant dispersion 
relation can directly be read from these expressions as 

Det \{ri^fik"kl^)T]^'' -k^" +X^'''(k)\ = (29) 

In this form it is explicitly formulated in terms of the undisturbed properly prescribed rela- 
tivistic plasma distribution function fso{p), holding in a homogeneous plasma, i.e. for wave- 
lengths substantially shorter than any inhomogeneity thus imposing restrictions on its appli- 
cability to long wavelength perturbations of the shock transition layer. 

The above theory has been formulated in terms of the 4-vector potential A^ which is 
subject to proper gauging. The most appropriate choice is the Lorentz-gauge k^^Afi = 0. In 
addition one may fix the direction of the wave vector k = Then k^ = (a)/c,k). 

This reduces the problem to the three independent polarisations: one longitudinal and two 
transversal polarisations with unit vectors e^,exi,fi±2- The explicit expressions for the po- 
larisation unit vectors are 

^/\r]apk^\e'^ = {k,,co/c,0,0), e^i = (0,0, 1,0), ex2 = (0,0,0, 1) (30) 

In this orthogonal system the 4x4 tensors D^' ,X^' reduce to 3x3 tensors Dij = e^^D^iveJ, 
Xij = A^v^J, and r\ij = erej, with i, j = 1 , 2, 3. As a result the dispersion relation becomes 
just a different version of Eq. ( |26[ l 

Det[D,y(A:)] = Det[r|f,yk^'k''^^iJ - kkj+Xijik)] = (31) 

In order to calculate the elements of Djj it is convenient to normalise the undisturbed dis- 
tribution function taking into account that the phase space volume is covariant: fso{p) — >■ 



^ _ ki^py+p^ky riafik^kf^P^Pv 



(riapk^pPy- 



(28) 



We use the flat metric tensor in the form r)f,v = T\^^ = diag(l, — 1, — 1, — 1). 
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fsQ (p) [Ns/y[^)]8 [p^ - mjC7(p)] , where = (m,c7(p) , p), p,- = p^" ■ e,-, and f,Q (p) is the or- 
dinary one-particle phase space distribution, which depends on the 3-momentum p, / dp/jo 
(p) = 1, and 7(p) = ^1 + Ipp/m^c^ is the internal gamma-factor. Then the reduced polar- 
isation tensor reads 



Aij 



-co: 



7 



dp/.o(p) 
7(P) 



riij + 



{k-k)pipj 



(32) 



(A more conventional form of this general expression can be found, e.g., in Eq. (8) of |Bret| 
|et al|2010| . Calculation of its components is straightforward but the result depends on the 
model. The non-zero components for a moving species s are (dropping the superscript s) 
Xi/, A^i XI, Ax2,±2, ^Li.e = ^( i_2- ^'^'^ counterstreaming cold beams and electromagnetic 
(transverse) modes they have been given explicitly (e.g. ,Silva et al,2002, |Achterberg and| 
|Wiersma|2007l l. 



Shock-relevant Wave Instabilities: The Fluid Picture 

For several reasons only the unstable solutions of the dispersion relation are of interest. 
They yield instability (growth, interaction, wave cascade, turbulence, and absorption) of 
plasma waves as being fundamental to shock physics. The astrophysical interest in all these 
plasma waves and instabilities which may generate them is rather limited; it is dictated 
by the observation of Cosmic Rays and radiation, with shocks themselves considered just 
mediators of both. Astrophysical interest in the very shock structure and formation (taken in 
this order) has been reluctantly arisen from the as well with reluctance accepted insight that 
probably neither particle acceleration nor its consequence, the generation of radiation, will 
be understood without at least some understanding of the principles of shock formation. 

Shock formation on its own is a consequence of instability. Firstly, shock waves result 
from instability of waves in collisionless plasma when these waves steepen by nonlinear pro- 
cesses, whether growing in the absence of dissipation or being pushed to grow by external 
drivers. Secondly, these waves are involved in the generation of dissipation, shock thermal- 
isation, and reflection and acceleration of particles with the latter two a collisionless way of 
dissipating large amounts of the upstream energy supply. Finally and, for relativistic shocks 
most important, a number of plasma waves seem to sign responsibility for amplification 
or even generation of the shock-required strong downstream magnetic fields by generating 
current flow in the shocked plasma environment of the shock, i.e. the apparently very high 
compression ratios Rs which the contemporary diffusive acceleration theories require, much 
higher than any gasdynamic and magneto-hydrodynamic theory can provide, even when 
pushed to its extremes. 

Because of this reasons, relativistic shock theories in the astrophysical context have in 
the recent years turned to investigate some shock-related linear instabilities. These are, in 
the first place, electromagnetic instabilities driven in the presence of the high energy par- 
ticle component. It is commonly assumed that the upstream flow is cold, an assumption 
that is good for at most mildly-relativistic temperatures T <C m,c^ and high Mach num- 
bers with F 3> 1 . In this case the distribution function degenerates to become a Dirac-delta 
function, and the fluid dynamical picture can to first approximation be used to investigate 
the relevant linearly unstable plasma waves. One distinguishes electrostatic (longitudinal, 
k||E, wave vector parallel to wave electric field) and electromagnetic (transverse, including 
mixed) modes. The last decade has seen high activity in analytical and numerical investiga- 
tion of these linearly unstable relativistic modes (sometimes also claims on their non-linear 
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evolution; as for an incomplete list one may refer to Achterberg 1983; Yoon and Davidson 
\19WP, 'Califano etal'1997'. 'Kazimura et ari998', "Zweibel 2002, |Bell|2004( |Blasi and Amato 
|2008i Bret et al 2010; Lemoine and Pelletier 2010 , 2011 ». 

Since the relativistic treatment of the full kinetic theory of wave growth is complex, the 
main analytical efforts have focussed on the approximate fluid picture which to some degree 
is equivalent to assuming a 'waterbag model' for the unstable set of particle distributions 
(as first used by lYoon and Davidson|1987] in an investigation of the non-relativistic Weibel- 
filamentation mode). There have been claims that the 'waterbag model' introduced spurious 
effects ( |Bret|2009] l but these are physically less important subtleties. The main problem with 
the fluid picture is that it applies only to a subclass of modes; it yields only 'reactive' and 
'non-resonant' instabilities and ignores thermal spreads of the resonances. 

In the cold non-magnetic (Oi = 0) beam approximation the beam distribution is Fh{p) <x 
5{px — r^inc^)S{p±) where ~ Yb is the beam (index b) Lorentz factor, and the beam 
propagates in jc-direction with px the transverse (cold beam) momentum p = {px,P±),Px = 



mcYhPx' ~ y ^ ~ ■ beam-fluid polarisation (susceptibility) tensor then becomes 

(33) 



ikPf + kjlif')c {co^-kV)l5!'l3f 

5ij + 



W - ckifif {CO - ckijij 



h\2 



For the cold-beam distribution the integral including the factor in front on the right-hand side 
then simply becomes a factor of — (Oft/o))^ before the expression in the brackets, where 0)^ 
is the beam plasma frequency. This susceptibility is to be added to the cold fluid suscepti- 
bility — {cOp/co)^, with co^, = tt)^(l +IJ.), and jU = nie/nii. The resulting dispersion relation 
has longitudinal electrostatic solutions, as Langmuir waves, and one short wavelength, small 
phase velocity co/k<^c electromagnetic very low frequency solution, the transverse Weibel 
('filamentation') mode for — s- 0,^, ^ with electric field along and magnetic field per- 
pendicular to the beam (and k^, therefore called filamentation). The growth rate is 



Im(ffl) = C0b/y\+ ajjklc^ -s- coi, for k-^ > cOp/c (34) 

Since the initial (upstream) magnetic field was zero, this instability generates a weak about 
zero-frequency magnetic field by filamenting of the microscopic beam currents (for the first 
astrophysical a pplications cf. Medvedev and Loeb,,1999, .Gruzinov and Waxman, 1999a, 
Gruzinov]|2001 Wiersma and Achterberg||2004[ Lyubarsky and Eichler||2006[ Achterberg] 



and Wiersma'2007'; 'Achte rberg et al|2007 1. This is a well-known process since the pioneer- 
ing work of Fried ( 1959). It has been elaborated in early analytical ( Davidson et al| 1972l l and 
numerical |Lee a nd Lam pe|1973|l work and was observed also in laboratory laser-plasma 
experiments |Tatarakis et al 2003 | l. Clearly, in the zero-temperature approximation, any 



temperature effects are excluded, and the original Weibel instability l |Weibel|1959^ caused 
by temperature anisotropy is not recovered here. Moreover, this instability is of very short 
wavelength and thus can generate magnetic fields only on a short distance of the order of 
the existence of the particle beam. Hence when the upstream beam penetrates the shock and 
is thermalised, the field is created just inside the thermalisation length, i.e. inside the shock 
ramp transition, and the already mentioned problem arises of how a magnetic field can be 
transferred to the far downstream. On the other hand, if the shock accelerated (cosmic ray) 
energetic particles in the downstream can be assumed to represent a beam then the magnetic 
field would exist all over their confinement scale. 

Among the longitudinal electrostatic instabilities, Langmuir modes are of little impor- 
tance compared with the electron Buneman mode ( jBuneman|1958[[T959[ ) which is the most 
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active inside the shock ramp current region to which it is restricted. The Buneman mode 
causes violent thermalisation of the electrons and has profound nonlinear effects resulting 
in the production of large numbers of electron holes with signatures in real space as well as 
in phase space subsequent electron heating in the flow direction and electron acceleration. 
The latter is of interest because the phase space holes accelerate a substantial fraction of 
electrons into a very cold electron beam while the bulk electron component is substantially 
heated up to temperatures of the electrostatic hole-potential drop, which is a large fraction 
of the initial energy stored in the current flow (cf., e.g., Newman et al 2001 j . 

The secondary cold electron beam escapes along the Weibel-magnetic field into the 
downstream plasma (note that it cannot escape far upstream because of the relativistic shock 
speed and the retarding shock potential which inhibits upstream acceleration of electrons) 
where it serves to generate secondary electron holes via the Buneman instability. The beam 
continuously regenerates itself in the holes while it heats the surrounding electrons along 
its path in the magnetised downstream matter up to a distance where its energy becomes 
insufficient to further feed the generation of electron holes. The Buneman instability (or its 
magnetised equivalent, the modified two-stream instability) is probably the main responsible 
for the generation of a hot electron gas surrounding the shock transition which becomes 
visible in the emitted high energy radiation spectrum. It might also be the main responsible 
in driving the Weibel mode via the electron temperature anisotropy it causes and which will 
grow when the filamentation instability ceases. This whole process is continuously fed by 
the relativistic or ultra-relativistic upstream flow. Thus the grossly unexplored combination 
of a Weibel-like instability and a Buneman-like instability seems to be of vital importance 
in the physics of relativistic shocks. 

This picture requires understanding how the magnetic field can be amplified or generated 
near the shock over distances much larger than the shock ramp transition. Since it is believed 
that the Weibel instability just generates short-scale magnetic fields in the region of the shock 
transition only, one is on the search for mechanisms capable of producing magnetic fields on 
scales larger than the electron or ion inertial lengths on which the Weibel instability grows.. 

Steps in this direction have been proposed by |Bell| ( [2004| l who assumed that the shock- 
accelerated Cosmic Rays signed responsible for the generation of return currents in the 
ambient flow. Once these undergo fluid instability they can non-resonantly grow and as 
well generate about zero-frequency magnetic fields. This mechanism has become subject to 
intense investigation (cf., e.g.,|Reville et al'2Q06', "SOPS'. 'Blasi and Amato'2008;'Amato and 



Blasi|2009j.Riquelme and Spitkovsky 2009 ; Ohira et al 2009 ; Bykov et al.201 lb. ; Sironi and 



Spitkovsky|20T T\ and others) both analytically and numerically. 



Electron holes have been predicted first by Bernstein et al 1957 as kinetic structures both in real and 
momentum space. Their spatial scales are between a few 10 and a few 100 Debye-lengths, identifying them 
as micro-scale features which move at speeds very close to the bulk current flow velocity. They cause non- 
gaussian organisation of the electron distribution in phase space (for textbook accounts see also, e.g., |DavId^ 
|son|1972l [Treumann and B aumjohann'1997j. The statistical theory based on coarse graining of the phase 
space has been given by Dupree 1 1972 1982). Electron holes resulting from current instability have been 
investigated analytically and numerically by (cf., e.g., 'Berk and RobertsU967"Scham el|1972|[i975[|Schamel| 
|and Bujarbarua 1983, Schamel and Luque 2005, Luque and Schamel 2005 and many others). Their most 
important properties are (i) that they contain strong localised electric potentials of the order of the available 
plasma-kinetic energy, (ii) reflect and accelerate electrons to energies much larger than thermal, and (iii) trap 
the low energy component in the hole potential well. More recent numerical simulations (for the most sophis- 
ticated examples see Newman et al 2001 2007 , Main et al 2006 Goldman et al 2007 1 have investigated their 
configuration space and phase space dynamics showing that electron holes are an inevitable consequence of 
plasmas which carry strong electric currents in which case the holes completely determine the micro-scale 
dynamics of the plasma. 
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In an ambient magnetic field B = B^z, in the ion-inertial (Hall) regime cOd <ti CO <^ cOce 
of unmagnetised ions and magnetised electrons, the Weibel-filamentation growth rate of the 
wave with wave vector perpendicular to both the beam velocity and magnetic field becomes 



Im(o)) = co^jcos 9B„/y 1 + co/jcos 6b„ for k, > Wp/c (35) 

In the limit of cos dyB = 0,k ^ k,, with n = kc/(0, the dispersion relation reads 

(e;_„2 + X,) (ei-n2 + ''AJ-e2 = 0, ei^l-4 + -4' ^2^— (36) 

and has as solutions the growing ordinary Weibel filamentation and a whistler mode. These 
waves grow on the electron inertial scale = c/cOg thus being of rather short wavelength. 
We will brief on the Weibel mechanism in the context of numerical simulations. 

|Lemoine and PelletierH20I0[|2011[ l discussed a number of other instabilities of which 
they assume that it will be possible to excite them in the shock-upstream region: whistlers, 
Alfven modes. Bell modes and their relative mutual limitations in dependence on the as- 
sumed ratio jii of electron-to-ion mass. From their investigation it seems that other waves 
grow faster than the filamentation instability and thus would outrun it being of greater im- 
portance (even though not generating magnetic fields implying that one of the big problems 
in ultra-relativistic shock physics remains unsolved). However, a wave dominates over oth- 
ers not necessarily if its growth rate is larger. Wave growth requires in addition (i) that the 
source of free energy exists long enough for the wave to grow to sufficient amplitude, (ii) 
that the growing wave starts from a sufficiently large thermal fluctuation level, and (iii) that 
the growing wave does not start decaying into other wave modes at an early stage. Neither 
of these conditions has so far been checked properly. For instance, in a cold plasma of zero 
temperature the thermal fluctuation level of electrostatic low frequency modes is just on the 
quantum level. On the other hand, the magnetic fluctuation level might be substantial such 
that the Weibel mode may grow to large amplitude already before the electrostatic modes 
have started. In non-relativistic reconnection it has recently been shown ( |Treumaim et al| 
|2010[ l that the Weibel instability selects the largest fluctuation level for growing. 

Whether the conclusion that the Weibel filamentation mode does not grow is correct or 
not cannot be decided easily. Our previous discussion showed that the upstream foreshock 
(precursor) region of the relativistic shock is extremely narrow and, therefore, any unstably 
excited waves will be spatially restricted to the vicinity of or very close to the shock ramp 
itself. Speaking about superluminal waves is relevant only for the phase speed which does 
not transport energy. Since the group speeds of the waves will always be less than c, the 
energy of the unstable waves will also be confined to the shock ramp and must accumulate 
there. This implies that the waves very rapidly reach their nonlinear state, and linear analysis 
breaks down. All kinds of wave-wave and wave-particle interactions take place after this 
happens, and it cannot be said that one or the other mode does not grow. The most probable 
case is that the accumulation of wave and particle energy at the shock ramp will provide 
a very rigid barrier for the upstream flow and cause rapid thermalisation inside the ramp 
with possible broadening of the shock ramp in the downstream direction such that the shock 
becomes a broad downstream transition region that will differ substantially from the narrow 
ramp obtained when investigating low-Mach number or non-relativistic shocks. 
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Shock-relevant Wave Instabilities: Kinetic Theory 

It has been recognised that a pure fluid approach to the unstable modes around relativistic 
shocks cannot properly account for the main problems. This has stimulated the investiga- 
tion of the kinetic theory of unstable waves with the emphasis on long wavelength waves. 
Because of the complexity of such investigations one should in principle turn to numerical 
simulations. In this section we briefly review the most recent analytical and numerical work 
only, restricting to the main most recent proposals ( |Rabinak et al||2010, ; Bret et al^,20r0) ). 
As for the case of the fluid theory these focus on waves in the upstream (foreshock or 'pre- 
cursor') region of relativistic shocks with the main interest in the generation of large-scale 
magnetic fields as the result of cosmic ray beam-upstream flow interaction (being of most 
interest in cosmic ray acceleration theory, cf.,e.g., Bell 2004, 2005[[Katz et al|2007[|Keshet| 
|et al||2009[ [Bykov et al|2011b[ and others), assuming that the Cosmic Rays have Lorentz 
factors 7 ^ F and can run far ahead of the ultra-relativistic shock. Hence, this theory is 
subject to the same restrictive comments as the above fluid approaches. 

The basic kinetic theory of relativistic beams passing through a cold plasma has been 
given long ago ( [Akhiezer et al|1975| l. Such streaming can excite magnetic waves both res- 
onantly and non-resonantly. The long-wavelengths resonant waves are Alfven waves which 
grow when the particle streaming velocity in weakly magnetised plasma exceeds the back- 
ground Alfven speed Vj, > Via- The growth rate in this case is well known (e.g., |Blandford| 
andEichler|1987!|: lm{co)res « COc.cr[{kvh / Co) - l]Nb/N, with N = I,-7V;. 



Some general conclusion on long wavelength purely magnetic Weibel modes can be 
drawn ( |Rabinak et al|[20T0l l including a drifting cosmic ray component with distribution 
function Fi{p) as function of the relativistic particle momentum p of species /. The total 
plasma frequency is given as the sum 

^P-E^NE"^, N>=jf/\F.{p) (37) 

including the integral over the relativistic distribution function F/ (which is a somewhat 
unusual definition of the partial density A',). It includes an average over the inverse (internal 
cosmic ray beam) Lorentz factor Yi(p) which results from the relativistic expression of the 
mass in the denominators of the partial plasma frequencies^ 

The cold plasma dispersion relation with background flow and cosmic ray beam in- 
cluded becomes a sixth-order polynomial. Only two of its solutions yielding instability (in- 
dices II , ± refer to the direction of the cosmic ray 'beam' v;,, and (• • • ) means averaging over 
Fi) 

Im(ffl) = +©4(1 + 0)2/^2 c2)-i + (fc||c/a)6)2(7-2)] J, fc||(7-2)<^^<a)p (38) 

When a hypothetical beam spread of angle 4 is taken into account, [Rabinak et al| ( |2010| l 
find that instability exists below a critical spreading angle Ad < A Qcrit for waves of fre- 
quency O) :^ jSA:|| c under the condition PAQ <C0i,/c0p, with j6 = v/^/c. This has been shown to 



This formula results from Eq. |32|. It assumes implicitly that the phase-space volume element dx^dp^ is 
a Lorentz invariant. Accounting for Lorentz contraction, the spatial volume element transforms fro m the co- 



Landau 



moving (primed) frame to the observer's frame as dx^ =dx'^ /y{p), /(p) = 1^1+ /nfic^ (cf., e.g. 
[and Lifshitz| 1 975] Section 10). The momentum-space element, on the other hand, can be shown to transform 
as dp' = j[p)dp" which cancels the extra proper (internal) Lorentz factor in the spatial volume element, thus 
yielding dx^dp^ = djc'^dp'^^. Indeed, the phase-space volume element turns out to be a Lorentz scalar. 
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hold for all axi-symmetric cosmic ray beams. [Rabinak et al] ( |2010[ l note that these long wave- 
length modes are not the fastest growing modes. Shorter wavelength Weibel-filamentation 
modes will grow faster, as has been known since |Weibel| ( |1959] l. 

In application to relativistic shocks, one has (Ob /(Op ~ (yh 1^)^^ > with Ub the fraction 
of post-shock energy carried by the cosmic ray beam, and s ^2 the approximate power law 
index of the cosmic ray spectrum. 'Far upstream' it is assumed that 2> F. Then, according 
to the above estimates, the long wavelength modes grow at rate Im(a)) ^ {cOb /cOp)k^c. Since 
the wave should grow on time scales shorter than the deflection time of cosmic ray particles 
Trfe/ ~ ym/eB by an angle ~ the ratio of growth to deflection time determines that 
waves can grow for transverse wave numbers (in shock frame quantities) 

2 

k^c/(Op>{ax/ab)^-r-\ for 7b{x)<r(alr\ (39) 



In electron-proton plasma (Ob/(Op ^ \/(^bl^{7/I") ' , and for growing long wavelength 



modes it is required (Lyubarsky and Eichler 2006 i that F > 10^^O!/,/0.1. Smaller Lorentz 



factors suppress the growth of long-wavelength waves. In addition the minimum shock- 
frame energy fraction of the cosmic ray beam as function of distance must be less than 

These quite general results for cold beams of a certain angular spread have been gen- 
eralised to a more complete plasma model including a thermal spread of the particle dis- 
tribution (cf. [Bret et al|2010[ where a table is given compiling the maximum growth rates 
of the two-stream, filamentation, and so-called 'oblique' instabilities in a cold beam-plasma 
system for different parameter choices). This can be achieved by referring to a relativis- 
tic thermal Maxwell- Jiittner distribution function Fmj{p} °<: exp{— [/(p) — jitPx]} instead 
of the waterbag cold plasma delta-distribution. This is an isotropic-temperature thermal- 
equilibrium distribution function which holds in the relativistic isothermal case[^It is nor- 
malised to one if using the factor of proportionality <§/4;ry^(p)A'2[<§/y(p)], where ^ = 
mc^/T, and K2{x) is the modified Bessel function of second order. Generally, relativistic 
flows with non-relativistic temperatures T < inc^ which are of interest in relativistic-shock 
physics, will exhibit negligible temperature effects. The parameter space of the two-stream 
and Weibel modes has been investigated for the Jiittner distribution in order to elucidate 
any residual temperature effects and to determine the transition between two-stream and 
Weibel/filamentation instabilities with the known result that at non-relativistic flows the sys- 
tem is dominated by two-stream modes. For relativistic and ultra-relativistic flows the known 
and expected result is that the Weibel/filamentation instability family takes over. Though, 
since many parameters are involved (beam density ratio, beam spread, temperatures, Lorentz 
factors, mass ratio) the system of solutions becomes complex and non-transparent in the 
thermal Jiittner case. Clearly, in the bulk moving frame one has F = 1 . Hence, temperature 
effects may, in this frame, affect the evolution of other instabilities which, when growing in 
the direction perpendicular to the bulk flow, will survive and might also influence the con- 
ditions at the shock transition. Effects of this kind have not yet been investigated properly 
(see, however, |Dieckmann|2009| l . 

In view of the application to particle acceleration and the generation of radiation it is 
of particular interest to estimate the maximum achievable (in the above sense of the general 
self-similar theory, i.e. the rms) amplitude of the stationary magnetic field that is generated 



For a discussion of its validity see Dunkel et al 2007 . 
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by the Weibel-filamentation instability (here we do not enter into a discussion of the electro- 
static modes even though these might be necessary in preparing the conditions for growth 
of the Weibel instability). Such an estimate requires the precise knowledge of the saturation 
(stabilisation) mechanism of the instability which by itself requires the solution of the full 
nonlinear dynamics of the wave-particle interaction. Since this is an impossible task, a par- 
ticular mechanism has to be assumed. Quasi-linear relaxation (i.e. the random-phase mean 
field approach) is invalid in this context because no depletion of the beam takes place. 

The usual assumption is that the magnetic wave saturates by trapping the beam particles 
in the magnetic mirrors of the wave. This is the simplest choice from a large variety of 
possibilities like wave scattering, cascade, wave decay and other parametric interactions 
including non-linear wave-particle interactions. It supposes that one single wave number k 
grows fastest and dominates the growth of all other waves, clearly contradicting the self- 
similar theory assumption. Sticking to this mechanism, the trapping (bounce) frequency 
in a sinusoidal wave-mirror in the wave frame is (O,,- = ^Jevi,b„/mYi, ,,.c, where b„ is the 
wave magnetic amplitude. Hence, the bounce frequency in the wave well increases with 
wave amplitude b^, beam velocity vi, and wave number k and decreases with trapped beam 
Lorentz factor Assuming that the wave stops growing once the bounce frequency Wn- ~ 
Im(a)) becomes comparable to the growth rate (cf., e.g., Silva et al|2002| ), which yields for 
the saturation amplitude of the filamentation instability 

C^^{Im[co«]}2(y,,,//3,,,,),^ (40) 

The average indicated by the angular brackets ( ■ • ■ ) has to be taken solely over the trapped 
(bouncing) particle distribution Fbjr, and the growth rate is given by any of the above ex- 
pressions for the filamentation instability, depending on the chosen model. An example of 
the Weibel saturation field based on the magnetic trapping saturation assumption is given 
in Figure [7] It must, however, be stressed again that the magnetic trapping is just one ap- 
parently obvious mechanism of saturation based on the assumption that no other nonlinear 
effects set on earlier and the wave grows linearly until trapping in the wave wells cuts it 
off. Whether this proposition is realistic is not known. Trapping requires quite large wave 
amplitudes, a fact that may suggest that long before trapping sets on other processes like 
cascading may have taken over and either wave saturation ends up with smaller rms ampli- 
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tude or the spectral energy density is distributed over a larger range in wave number and 
frequency space. 

CR Modified Shocks: Bell-like Instabilities 

One particular subclass of electromagnetic instabilities is worth to be singled out. These 
are the Bell-like instabilities (originally anticipated though not worked out) by |Achterberg| 
( |1983[ l. They were recognised in their importance and proposed as a fundamental mecha- 
nism for generation of magnetic fields near shocks in the presence of Cosmic Rays only 
more recently by Bell (2004, 2005). They differ from the Weibel-filamentation mechanism 
which is based on beam-beam interaction (or temperature anisotropy). Bell-like instabilities 
are non-resonant electromagnetic very low frequency instabilities which rely on the pres- 
ence of cosmic-ray 'return currents' flowing in the shock background plasma environment 
and closing the presumable current that is believed to be carried by the shock-accelerated 
high-energy Cosmic Rays into upstream. Of course, they are in principle already contained 
in the general theory given above and have more recently been investigated based on general 
kinetic theory (Reville et al 2006| [20071 |2008l [Amato et al|[2008l [Blasi and Amato 



2008^ 



Amato an d Blasi||2009 1 and also with particle-in-cell (PIC) numerical simulations Ohira| 
et al|2009| ). 

The original proposal by |Bell| ( |2004[ ), following numerical simulation work by |Lucek| 
[and Bell| ( |2000[ l, was based on a fluid (MHD) model including a hypothetical cosmic ray cur- 
rent density Jcr = eNc'^h (with vi, the cosmic ray beam velocity) which the background 
plasma compensates (on its frame of reference) by an electron return current of Jret = 
—eNcr^, assuming that for reasons of charge neutrality this current transports the same 
charge density as the Cosmic Rays, and its velocity is the background plasma velocity in- 
duced by the Cosmic Rays and yielding the additional Lorentz-force term — {Jcr — J ret ) x B 
in the MHD-momentum equation of the upstream magnetised background flow in its proper 
frame. |Bell^ ( 200 4 ), working in the upstream frame, assumed a parallel shock such that the 
original upstream magnetic field is parallel to the shock normal. /|| is the free input param- 
eter generating a first order perpendicular magnetic field component Bx and perpendicular 
current Jx = a(y||/S||)Bx. The quantity a is calculated to first order as function of A = 
{^''c,cr) the inverse product of wave number and upstream cosmic ray gyro-radius. This 
yields a linear upstream frame dispersion relation O)^ — k^V^^ — R[l —a{co,?i) — a/kvi,] = 
with R = /miN which is a cosmic ray modified upstream dispersion relation contain- 
ing Alfven waves. For large cosmic ray velocities v;, the last term in the brackets can be 
dropped, in particular as the wave frequency is small, CO ~ 0. The Cosmic Rays turn out to 
drive a magnetic wave field in the upstream plasma resulting from a non-resonant instability. 
For a power law cosmic ray momentum distribution F{p) cx p^'^ in the interval pi < p < P2 
the quantity a can be calculated. Bell,(2004) derived the dispersion relation 

CO^ - fcVii ± Cvi(l - ai)k/r,^,n = 0, where ^ = rc.cri \BJ\i l/vlnnN (41) 

and ai — 1 for Ai = (A:rc,„i)^' ^ 1 while ai — )• for Ai <C 1. For Ai = 1 a stable Alfven 
wave is obtained in the uninteresting case ^v^ <C V^^ , when the Cosmic Rays do not carry 
currents. In the opposite case the dispersion relation is (O^ w ±^(V|(,/4r(..cri)^(l — ^Tti) and 
instability arises. These are purely current driven non-resonant non-firehose modes, very low 
frequency magnetic oscillations driven by the cosmic-ray current of particles which on the 
unstable short scales are non-magnetised. Conditions under which this theory is applicable 
have been discussed in depth by jBykov et al| ( (201 la[ l. 
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The most relevant instability is found for moderately short wave lengths 1 < kr^cri < 
(vfj/V/^iYC,. The maximum growth rate becomes Im(a))|m w j^^MOci.cn where = 
vi,/Va\, with maximum unstable wave number w {vb/yAx)^ /'^fccri, and the predicted 
saturation wave amplitude 5B/B\ <x ^/^i,APcr / miNiV^/^ w 300 is quite large. Note that the 
assumption is that the initial upstream field is not zero. Nonetheless, in a very weak ini- 
tial upstream field Bi the saturation level may still be quite low. In favour of the Bell in- 
stability one should, however, note that the saturation level is only weakly dependent on 
the background field; it is determined by t he cosmic ray pressu re at the shock surface: 



SB-^/Sk ~ {Vsh/2c)Pcrj. On the other hand, Reville et al (2006 1 have shown that plasma 



temperature effects lead to further reduction of the growth rate and also affect the most un- 
stable wavelengths. Still, Cosmic Ray diffusion seems to be amplified substantially, coming 
close to the Bohm limit, as MHD simulations by the same authors suggest (Reville et al 
2008V Moreover, numerical PIC simulations performed recently fRiquelme and Spitkovsky 



2009 ) seem to indicate that due to early nonlinear saturation of the instability the amplifi- 
cation which can be achieved is much less than Bell's above estimate, barely exceeding a 
factor of 10. Related numerical simulations by [Niemiec et al| ( (2008j ) with a rather different 
setting had already previously suggested even lower saturation levels, barely exceeding the 
background magnetic field value. 

There are additional reservations concerning the Bell mechanism with return current in 
the upstream region. The validity of the magnetohydrodynamic approach in the Bell model 
(where the return current is carried by background electrons that move together with the 
Cosmic Rays, which is possible only when the shocks are parallel) has nevertheless been 
confirmed by kinetic calcula tions (|Reville et al||2006[ [Amato and Blasi||2009[ l and by the 
PIC simulation of |Ohira et al| l |2009| l. The flow of the currents in the upstream region implies 
that the (highly relativistic) shock approaches with light velocity such that it remains ques- 
tionable whether the comparably slow waves will have time to grow fast enough to reach 
any substantial amplitudes for saturation before being overcome and digested by the shock 
ramp. 

Also, the maximum growing wave numbers remain to be large since the cosmic ray 
Mach number S> 1 is large and, thus, the unstable waves are of quite short wave- 
lengths. This was believed to be in contrast to the resonant streaming instability modes (cf., 
e.g., Lagage and Cesarsky'1983'); however, [Amato and Blasij ( |2009 ) have shown that both 
types of waves have maximum growing wave numbers krc^cr ~ 1 and comparable though 
small growth rates at long wavelengths (see Figure [8|; [Reville et al| ( |2006] l performed sim- 
ilar calculations for the relativistic kinetic case including temperature effects which yield 
essentially a similar result thus confirming that growth will occur for wavelength below the 
gyro-radius of the resonant cosmic ray particles though closer to the gyro-radius than in the 
non-relativistic case. This result was also confirmed in the simulations of |Ohira et al| ( |2009l ). 
One would not only like to have long wavelengths but also a theory for mostly perpendicular 
shocks because in fast streams with weak magnetic fields, i.e. Oi <^ I the shocks are very 
close to being perpendicular. 

[Amato and Blasi| ( |2009[ l have checked the competition between the non-resonant (Bell- 
like) and resonant modes in a fully kinetic approach for the Bell settings and for the [Achter^ 
[berg ( 1983 1 case when the cosmic ray current generates a compensatory drift in the upstream 
background, i.e. for the density Ne = A',- +Ncr and zero current NjVi = N^v^ conditions. In the 
latter case the dispersion relation agrees with Bell's only for Nc- <C A^;. They extended their 
calculation to a cosmic ray power law spectrum <x p^^ with 4 < j < 5. It turns out that the 
relative importance of the modes is determined by the parameter Z = ^(v^ /Va)^. For Z 2> 1 
the non-resonant Bell-like mode grows almost non-oscillatory and at a fast rate, dominating 
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Fig. 8 Comparison between the kinetic reso- 
nant and non-resonant cosmic ray return cur- 
rent driven magnetic upstream modes for the 
case when the non-resonant modes is the more 
important zero frequency oscillation, i.e. in 
the parameter range Z S> 1 (after |Amato and| 
|Blasii2009j l. The numbers chosen are Vi = 10* 
kms ', Bo = IflG, Ni = 1 cm"^. The calcu- 
lation assumes a cosmic ray acceleration effi- 
ciency of J] = Pcr/injNiV^ = 10"' and a max- 
imum cosmic ray momentum Pniax ~ 
Frequency (O and growth rate Im(B are given 
as functions of wave number kr^cr- Long 
wave length grow at similar rates for both 
modes while the non-resonant Bell-like mode 
dominates for short wavelengths. 



the resonant mode. Conversely, for Z ^ 1 the resonant modes dominate. Figure [8] shows 
a comparison of the two modes, resonant and non-resonant for the case when Z 1. The 
growth rates of resonant waves maximise at krc^cr ~ 1 while the non-resonant growth in- 
creases deep into the short wavelength domain, and in agreement with Bell's prediction the 
non-resonant mode takes over at short wavelengths kr^cr > 1 and the wave becomes about 
zero frequency. For long wavelengths krc^cr < 1 there is no difference between the growth 
rates and frequencies of the two modes. 

The Oblique Case 

It has already been noted that Bell's case, even though supported by kinetic theory, is 
marginal as it does not include any obliquity neither for the magnetic field nor for wave prop- 
agation. It holds just for a strictly parallel shock while it is suggestive that ultra-relativistic 
shocks will rather be either magnetically oblique or perpendicular. 

Obliquely propagating waves should be important in shock physics and shock-particle 
acceleration. This is known from non-relativistic shock theory and simulation (cf., e.g., the 
comprehensive account of the various plasma waves and their importance in all aspects of 
nonrelativistic shocks given in |Balogh and Treumann|201 l| l. In relativistic shock theory the 
importance of obliquity of waves has been noted as well (see the wordy but non-transparent 
discussion in Bret 2009) . However, in view of the very large number of sources of free en- 
ergy in the interaction of relativistic streaming inhomogeneous magnetised/non-magnetised, 
current-carrying/current-compensated plasmas with their environment - as in the case of ex- 
ternal shocks - with inclusion of high energy accelerated particle populations, any simple 
listing of the possible instabilities of all kinds must necessarily be incomplete and erratic. It 
just contributes to confusion. 

In non-magnetised flows (upstream) obliquity refers to the angle between upstream flow 
and wave vector k. In the presence of a background magnetic field however weak or once 
a magnetic field B is generated (as, e.g. in the cases of Weibel's and Bell's instabilities) it 
refers to the angle between B and k. Growing magnetic modes in this case should not have 
purely magnetic properties like non-resonant Bell modes; because of the implied obliquely 
directed electric field component, obliqueness rather introduces density fluctuations, for 
large- amplitude Alfven or the wanted Bell modes (at least in second order) and, thus, the 
waves become compressible, belonging to the family of magnetosonic waves. 
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Fig. 9 Growth rates of the parallel Bell mode 
(same as in FigurejsJ compared to the growth 
rate of the long-wavelength mode that is ex- 
cited in the presence of short wavelength 
Bell-mode turbulence. The calculation uses 
the following initial parameters: Bell-mode 
wavenumber of maximum growth kor^o = 
100, cosmic ray spectral power law index a = 
4, mixing length paramter = 5, and mean 
free path parameter 7] = 10 (after |Bykov et al| 
|201 Ib^ . Shown is the mode with the largest 
growth rate. Clearly the turbulence driven 
long wavelength mode has increasingly larger 
growth rate than the Bell mode in the long 
wavelength regime. 

10-3 10-2 lO-l 10° 10^ 10^ 



In initially even weakly magnetised upstream flows the main currents carried by the 
energetic particle component in their gyration motion. These currents flow perpendicular 
to the magnetic field; they cannot be compensated in the manner proposed in Bell's orig- 
inal mechanism ( ^Bell 2004^) because the magnetised electrons do not follow the gyrating 
energetic nucleons. They close via field-aligned currents carried by electrons which are ac- 
celerated along the magnetic field out of the inflow-plasma background. While the gyrating 
nucleonic-cosmic ray-current rings excite both electrostatic modified-two stream instabil- 
ities and, through the resonant ion-cyclotron maser mechanism, upstream low-frequency 
electromagnetic cyclotron instabilities (cf., e.g., |Amato a nd Arons'2006), the ionic equiva- 
lent to the celebrated resonant electron-cyclotron (cf., e.g., Treumann 2006) maser, the field- 
aligned electron closure currents are distributed over the entire cosmic ray cyclotron orbit. 
They can undergo other instabilities both electrostatic and electromagnetic. [Bell| ( |2005) gave 
a first inclusion of transverse currents, and |Bykov et al| ( [2009[|20H"b| l included a magnetised 
cosmic ray current, showing that short wavelength perturbations affected the longer wave- 
length instability. First numerical PIC simulations performed showed that the amplification 
of magnetic field in the parallel current Bell instability is lower than analytically estimated 
( [Riquelme and Spitkovsky|2010 |. Under the assumption that some cosmic ray particles re- 
main magnetised and drive a perpendicular current flowing closer to the shock they however 
suggest an additional short-wavelength amplification of the Bell-like magnetic field. This 
will be discussed below in the simulation section. 

Transition to Long Wavelengths 

The point which interests us here is how Bell's mechanism extends into the long-wavelength 
domain obtaining amplified magnetic fields on the larger scale, as required by the diffusive 
shock acceleration mechanism. A step into that direction has been gone by [Bykov et al| 
pOIIa|b[ l. These authors propose that the short scale magnetic field perturbations generated 
by the Bell instability provide a turbulent background for large-scale mean fields to grow. 
The reason is that the current carried by the cosmic ray component responds to the presence 
of the short scale turbulence, shorter than the cosmic ray gyroradius which, in the quasilinear 
(mean field) approach, scatters the Cosmic Rays diffusively and so diffuses the cosmic ray 
current. In the last consequence, this leads to damping of the short-scale kr^cr > 1 modes and 
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growth of long-scale kr^cr < 1 modes with propagation parallel to the initial magnetic field 
on the expense of the short scale modes. This dominance of growth of the new (diffusively 
excited) long-scale modes over the Bell modes is shown in Figure|9] 

The response in the mean field equations is modelled by assuming diffusive scatter- 
ing of the current carrying Cosmic Rays at the short Bell fluctuations which impose a fi- 
nite (anomalous) mean free path on the Cosmic Rays. The latter heuristically modelled as 
iinfp = TlfcO, being proportional to Kcq, the initial gyroradius of the Cosmic Rays in the undis- 
turbed field Bp, assuming 7] > 1. T he extreme Bohm limit (fastest diffusion) would corre- 



spond to r] = 1. Bykov et al ( 201 lb i use the Bell velocity fluctuation spectrum 4;rp |v(k| = 
|b(k)pA:i/|A:||| as function of the magnetic fluctuation spectrum |b(k)p in Bell's modes. 
This fluctuation spectrum is understood as a given background turbulence. Averaging the 
dynamic equations and the induction equation properly in its presence produces a prefac- 
tor (1 + (Cf) on the mean Lorentz ('ponderomotive') force term in the dynamic equations, 
where K, is a turbulent coefficient that is responsible for the transfer of energy into the long 
wavelength regime. It can be expressed through the average Bell-mode magnetic fluctuation 
spectral density (b^) and Tcoir ~ '^^^/ko\/ the correlation time, as 

I 

'^corr f . . ( 



with z along the undisturbed magnetic field Bq, h) = AnJcr/cBQ the wave number of max- 
imum growth of Bell modes, mixing length parameter ~ 5 for short wavelength Bell 
modes, and k\ the (mean) threshold wavenumber 1 < kvcQ < ki rco of Bell's modes, defined 
by ki(B) = {47z/c)(Jcr) where averaging is over the short modes, and /„• ~ ecNc-Ocr, with 
Gcr > B^/4K£cr referring to Cosmic Rays only. 

[Bykov et al| ( |2011b[ l obtain a general dispersion relation which in the 'intermediate' 
regime, where the wavelength of the unstable modes is longer than that of Bell modes but 
still shorter than the mean free path, i.e. r] ^ ' < krc^cr < 1 > yields a parallel propagating wave 
with right-hand polarisation and growth rate 

hi„, = bi„,{x + iy), lm{co)i„,^47ikVA{^{byBl))'^ (43) 

The polarisation of this mode is opposite to that of the left-hand polarised Bell modes 
(and thus offers its possible capacity of partially balancing the magnetic helicity of the 
generated magnetic fields). For even longer wavelengths in the 'hydrodynamic' regime 
krc,cr < 7]^' when the mean free path is shorter than the unstable mode wavelength, the sys- 
tem becomes diffusively coUisional and can be described by appropriately defined magneto- 
hydrodynamic equations. Their linear unstable solution s are both, left and right-hand po- 
larised, waves which grow at the same rate 



Im(ffl),,,rf K l^^^y kVA, ^ = cos > cos e,na, ~ - (44) 
Z TjK Bq J k Tj 



These very long wavelength waves are definitely oblique, and the maximum growing modes, 
for large T] , propagate about perpendicular to the magnetic field. 

The critical dependence of the unstably excited long wavelength magnetic fields on the 
finiteness of the mean free path ?i,„fp through the parameter 77 implies that field energy is 
dissipated. This dissipation goes on the cost of the driver, i.e. the Cosmic Rays via Bell's 
instability, but the dependence on 77 implies that some part of the cosmic ray energy is also 
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transferred to the background plasma, predominantly the ions in this case. [Bykov etal| ( (2008) ) 
have reviewed ion heating and heat transfer to electrons in shocks for the case when it can 
be modelled as kind of Coulomb collisions (justified approximately by the assumption of a 
finite mean free path Xmfp)- 

The work of [Bykov et al| ( [2011b| l ist the first quasi-linear (in the literal meaning of the 
word) attempt of describing the self-consistent generation of low frequency quasi-stationary 
magnetic fields in the shock environment - actually just in the upstream pre-shock region. 
It shows that a gyrating cosmic ray gas that by the Bell mechanism excites short-scale mag- 
netic fields in the lowest non-linear (i.e. quasi-linear or mean field) approximation is diffused 
by its self-generated short-scale wave spectrum, and the response of the upstream plasma to 
this diffusive partial depletion of the cosmic ray beam is to excite another branch of wave 
spectrum the scales of which are much longer than those of Bell modes. These waves es- 
cape scattering and can grow on the expense of the cosmic ray and Bell modes.This theory 
is interesting and promising in providing (still short though somewhat) larger-scale waves 
than obtained by simple linear Bell theory. It is a first and obviously efficient step toward a 
theory of inverse wave cascading from short into long scales as known in some models of 
turbulence where, as in this model, the wave energy-injection source is at the short scales, 
the scales of Bell's modes shorter than the cosmic ray gyro-radius. 



4 Simulation studies 



Most progress in genuine coUisionless shock theory is attributed to recent numerical simu- 
lation studies based on full particle, either PIC or Vlasov, codes. This is probably the only 
way of self-consistently treating the arising highly nonlinear problems when dealing with 
high-Mach number coUisionless shocks to which fluid theory is an insufficient approxima- 
tion. It should, however, be stressed that the astrophysical interest is by no means in the 
physics of shock formation and shock structure; its focus is on the capability of shocks to 
accelerate particles - on the one hand for reproducing the observed (or inferred) energy 
spectrum of Cosmic Rays, on the other hand for understanding the observed energetic pho- 
tons which are generated in those systems and which are believed to necessarily result from 
particles which have been accelerated at the shock. In this view, the relativistic shock struc- 
ture itself is just a mediator. The model underlying these expectations (cf., e.g., [Drury|1983[ 
[Blandford and Eichler|1987| [MaUcov and Drury|2001| l is the diffusive shock particle accel- 
eration model (DSA) which, for being effective, requires strong magnetic fields in the shock 
environment. We noted already that this implies that the shock self-consistently generates 
magnetic fields. Since the need for field amplification was realised, the main interest in rel- 
ativistic shock simulations focussed on magnetic field generation and amplification up to 
magnetic field strengths which substantially increase the classical MHD compression ratio. 
Following the initial work by Yoon and Davidsoii|i 1987^, this happened roughly one decade 



ago (|Nishikawa et al|1997[[Califano et al|1997[ 7998a|b[[kazimura et al 



1998; Medvedev ami 



[Loeb|1999[ [Gruzinov and W axman' 1999a ; Bykov and Uvarov 1 999 ; Nishikawa et al|2003) . 
These authors concentrated on the effects of the Weibel-filamentation instability. Their work 
caused a sudden inflation of activity in theory and simulation on the Weibel-filamentation 
mode which has not yet ceased at the time of writing (cf., e.g., [Dieckmami et al[2007[[Dieck-[ 
[mann|2009 ). It was later joined by Bell's instability, which however requires the presence 
of already otherwise produced Cosmic Rays. In the following we review recent numerical 
work in view of how relativistic shocks form and how they contribute to the needs of as- 
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trophysics - with the exception of high energy laser fusion plasmas, the only place in the 
Universe where relativistic shocks exist. 

Limitations 

Before proceeding, we note that all numerical simulations encounter a number of serious 
limitations. These should be kept in mind prior to naively applying the simulations to a 
given astrophysical situation. Published simulations are frequently not transparent enough 
such that their astrophysical validity can from time to time hardly be inferred. 

The first of these limitations is related to the number of particles that can be accessed in a 
simulation. This number is limited by the capability of following the paths of all particles in 
a given simulation box. Available computing powers allow treating only a comparably small 
number of particles of say (in the very best case) < 10'^ in total. This number determines the 
physical size of the simulation box, cells, and resolution. It can be compared to the number 
Nd = ^NX^ in the Debye sphere which must be large. No 2> 1, in order for the plasma 
picture to be valid. 

In relativistic shock plasmas of 7 ~ 10 keV <C m^c^ and A' ~ 10^^ m^^ for instance, 
one has A/j ~ 5 x lO-' m. No ~ 10'^ particles per Debye sphere, implying that not even the 
number of particles in one single Debye-volume could be simulated by present techniques. 
Thus, each 'simulation particle' is a macro-particle which represents a huge number of real 
particles; this supposes that some microscopic mechanisms are implicit to the simulation 
which cause the particles to be correlated (or clumped) together on the micro-scale for the 
entire simulation time Tj,„, such that they obey a common dynamic. One may assume that 
these are just the particles in the Debye sphere which are held together by electrostatic 
fields. Then the number of macro-particles contained in the simulation box is in fact the 
number of Debye-spheres which is also suggested by the requirement of charge neutrality. 
For sufficient statistics each simulation cell of linear size dc^us should contain a substantial 
number of such macro-particles; in the best PIC simulations these are about ~ 100 macro- 
particles per cell corresponding to a linear cell size of dceiis ~ 5 x 10^ m. Managing the 
calculation of ~ lO'^ macro-particle paths then allows designing a linear ID-simulation of 
Nceiis ^ lO'" cells corresponding to a physical scale of Lid ^ 5 x lO'^ m ~ 0.2 pc. In higher 
dimension these numbers (10 particles per cell in 2D and 3 in 3D) for reasons of reliable 
statistics require a reduction of cell number (to 10^ in 2D, and lO-' in 3D) and thus apply to 
shorter linear physical scales. 

Relativistic shocks move at high speed V^i, ~ c across the box. They quickly encounter 
the boundaries where the results become affected by boundary effects. This restricts the as- 
trophysically useful simulation time. On the expense of the box size, the boundary problem 
can be softened by producing the shock in a collision of two counterstreaming plasmas in- 
stead of a reflecting wall and leading to the generation of two shocks; in the case of external 
shocks, however, this method suffers from the necessity of considering two media of differ- 
ent physical properties: the stream and the environment. Another possibility is filling the box 
with a medium of artificially high dielectric constant e I. This, on the other hand, affects 
the growth of instabilities and wave propagation, thus causing other unwanted effects. In 
addition, | Vay | ( (2008) recently noted that the standard method in numerical simulations used 
to advance the particle momentum in time, the so-called Boris pusher, may break down at 
relativistic particle velocities because the split of the electric and magnetic components of 
the acceleration inherent to it causes a spurious force on the relativistic particles. Though 
the effect of this spurious force is not easily quantifiable in highly variable (turbulent) fields, 
the interpretation of particle motion and the ultra-relativistic part of the distribution function 
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Fig. 10 Relativistic perpendicular shock for- 
mation in ID-simulation of a cold charge neu- 
tralised pair plasma of density = N- and 
upstream Lorentz factor F = 40. The shock 
is produced in the interaction between the in- 
flow beam (coming in from left flowing in +.v 
direction) and the beam that is reflected from 
the impermeable boundary at the right and 
flows in —X direction. Shown is the 4-velocity 
«, (top), which initially is ii F, magnetic 
field B^, and electric field £,.. The fields are 
normalised to Bq. The spatial coordinate is 
in terms of the initial gyro-radius r^o. The 
shock propagates to the left. It causes elec- 
tromagnetic precursor waves in the field com- 
ponents which slightly distort the incoming 
beam in front of the shock (the large turnover 
in H j to negative values) and subsequent beam 
plasma thermalisation with downstream tur- 
bulence and production of a small number of 
energetic particles only up to twice the up- 
stream 4-velocity (i.e. up to F ~ 80). This ac- 
celeration is not seen upstream; it occurs sev- 
eral gyro-radii downstream. Compression ra- 
tio is « 2 (data taken from |Amato and Arons| 
p006) . 



require some care. This effect becomes crucial in strong electromagnetic fields; in the rather 
weak magnetic fields related to relativistic shocks it might be less important. 

Finally, the ion-to-electron mass ratio jU,e = miYi/nieYe causes troubles as it sets the 
scales of electrons and ions far apart. Since relativistic electrons must be assumed to be 
cooler than ions (because of their high radiation efficiency), the presence of the internal 
Lorentz factors in the above ratio make the problem even worse. For this reason, simula- 
tions in the relativistic domain are usually performed either for pair plasmas (and thus for 
mass ratio jU = 1), or using small mass ratios jU.,,-,,, <C of the order of at most a few times 
100. Simulations of non-relativistic shocks (cf., the review of shock theory in |Balogh and| 
|Treumann||201 1| | unambiguously demonstrated that the behaviour of collisionless shocks 
changes substantially when the mass ratio becomes realistic and that low-mass ratio sim- 
ulations can be quite misleading on what concerns shock formation and shock structure. 
Whether this also holds for relativistic shocks is not known. One may suspect that similar 
changes occur when it will become possible to increase the mass ratio also in relativistic 
simulations. 



4.1 Pair Plasmas 

The simplest collisionless particle simulation models use electron-positron pairs ±e in ID 
spatially. A typical yet highly idealised simulation of this kind l |Amato and Arons|2006[ l is 
shown in Figure[T0]for a strictly perpendicular relativistic shock of moderate upstream bulk 
Lorentz factor F = 40 and a = 1 (i.e., o = CT+ = 2), with a box of Nceiis = 1024 cells, 
dceiis = 0.1 and time step At = dceiis/c The simulation uses a reflecting impermeable 
wall, and the shock is caused by the interaction of the two counter streaming (cold) plasmas 




Collisionless Relativistic Shocks 



39 



in ±jc directions. As expected, one observes an electromagnetic precursor wave running 
ahead of the shock at the speed of light. 

The shock occurs where the incoming beam velocity is turned over into upstream di- 
rection and is accompanied by strong field oscillation and magnetic overshoot. Further 
downstream the plasma becomes thermalised with turbulent field and compression ratio 
w 2. Interestingly, a small number of energetic particles is generated downstream reach- 
ing a Lorentz factor F w 80. These particles do not reach the upstream region, at least not 
in the simulation time as their upstream directed velocity is not high enough to overcome 
the shock, and there is not the slightest sign of upstream particle acceleration. [Amato and] 
|Arons|p006 t argue that these are simply thermalized particles, estimating that the effective 
downstream temperature reached is T2 « \5mc^, which implies that the downstream heating 
is strong, and converts the downstream region into a region with relativistic temperatures. 
This, in turn, would imply that processes like particle creation/annihilation would ultimately 
come into play. It is nevertheless interesting that the tail of the downstream 'quasi-thermal' 
distribution is long enough to show the presence of quite energetic particles which have been 
accelerated in the downstream turbulence and not at the shock proper. 

Simulations in Higtier Dimensions 

The first higher dimensional PIC simulations with astrophysical applications have been 
performed by |Silva et al| p003[ l. These were 3D fully kinetic electromagnetic relativis- 
tic simulations of shock waves produced in the collision of two interpenetrating shells of 
electron-positron pair plasmas. The total number of particles whose paths could be followed 
amounted to just w 10**. Since the focus was on the generation of magnetic fields via the 
Weibel/filamentation mode, the shells were assumed to be initially non-magnetised with 
CT = 0, and the simulation was performed on a box of 25.6^ x lOA^ cells for a simulation 
time of TCOe = 150 and periodic boundary conditions. Here, as before, ?ie = c/cOe and cOe 
are the electron-skin length (electron inertial length) and electron plasma frequency, respec- 
tively. The inital Lorentz factors used were F w 1.17 (weakly relativistic) and 10.05 with 
thermal spreads 4/3//, w 0.085 and 0.01. Under these conditions the Weibel-filamentation 
mode grows, the current decayed into narrow filaments, and in the interaction region (which 
was on the size of a few /l^ ) a quasi-stationary magnetic field evolved in both cases, in the 
relativistic case within times w 150)7', afterwards saturating on an about constant level. 
These qualitative observations confirm the Weibel linear and nonlinear theories of the pre- 
vious sections, but quantitative application to astrophysics was hardly possible. 

A similar 3D-model of counter-streaming pair-plasma shells of equal density and bulk 
momentum being initially charge and current free has been used by [Jaroschek et al| ( |2004[ 
|2005| l again with emphasis on the Weibel-filamentation mode using a similar but slightly 
bigger grid with 5 x 10^ particles. The simulations use Lorentz factors F = 10 for inter- 
nal shocks and F = 100, the latter corresponding to the ultra-relativistic case, and thermal 
spreads jS//, = 10^^ and 10^*', respectively. In the ultra-relativistic case, which models the 
external shock scenario relevant in GRB afterglows, simulation times were Tfflg = 260, cor- 
responding to roughly 5 x 10^ time steps. In the mixing region of the two plasma shells the 
shock structure evolves into current filaments surrounded by toroidal magnetic fields. Grad- 
ually, neighbouring parallel filaments merge until saturation when the filaments are large 
and at sufficient distance separated by their proper toroidal fields. 

The instability is an oblique combination between the Two-stream and Weibel filamenta- 
tion instabilities but, for increasing Lorentz factor, evolves into a quasi-perpendicular Weibel 
mode with propagation perpendicular to the current-magnetic field structure, implying that 
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Fig. 11 Average magnetic (cj) and electric 
(e) energy density ratios in the collision 
of two pair shells for two different Lorentz 
factors F = 10 (red curves) and F = 100 
(blue curves) as function of time. Due to 
the stronger coupling of the two-stream and 
Weibel modes the weakly relativistic case sat- 
urates earlier and at a lower level than the 
ultra-relativistic case. The Weibel perpendicu- 
lar magnetic level is considerably higher in the 
latter case. Moreover, in saturation the trans- 
verse magnetic field energy density is five or- 
ders of magnitude higher than the longitudinal 
(data taken from |Jaroschek et al|2004) . 

100 ^ 200 
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the system becomes quasi-2D though some obliquity is necessary to keep the global system 
current-free. It is important to note that strict two-dimensionality would require the closure 
of currents somewhere 'outside' along the magnetic field and would thus cause different 
effects related to field-aligned current flow. 

Figure [TT] shows the time-evolution of the normalised energy densities in parallel and 
perpendicular directions to the counter-streaming shells. Here the volume-averaged energy 
is defined as (ct) = V^^ J dr^(S^/2jUo)/4(F — \)NmeC^. The weakly relativistic case (red 
curves) saturates early and at a lower level than the ultra-relativistic case. These fields are 
confined to the interaction region which is narrow, and the simulations do also suffer from 
the small boxes and the limited simulation time. Saturation is due to particle trapping and 
thus is diffusion limited. This yields life-times of trapped particles and magnetic field decay 
ofTdiff^ 1 0^ 0)7 a comparably short though reasonable time of ~ 10^ s, i.e. of the order of 
one or few days in real times. This time is not far from the afterglow time. One may notice 
that nucleonic components will shift the action to the ions thus extending all times by a 
factor which is a function of the mass ratio jj. = mi/m^. The electrons, after initially having 
contributed to seed Weibel fields become readily magnetised and behave passively while 
the ions will take over letting the electrons radiate in a magnetic field which will further 
increase, this time now on the ion time scale. 

Shock Structure 

Systematic investigations of the shock structure have been initiated by [Nishikawa et al| 
([2005 ) and S pitkovsky (20 05[l who performed 3D PIC simulations on perpendicular shocks 
in pair plasma. [Nishikawa et al| p005[ l used a comparably small simulation box with an 
85^ X 160 grid, = 4.Sdceii, and performed several independent simulations with a moder- 
ate total particle number varying between 5x10^ and 18 x 10^ particles (a number which 
distributed over the available simulation cells this corresponded to 27 particles per cell per 
species); they also applied periodic boundary conditions. The main difference their and other 
simulations as, for instance, those of |Spitkovsky| ( [2005] l which we are going to discuss in 
more depth, is that they separately injected fast spatially limited jets of different width into 
a given plasma set up at rest. They used both very narrow and wide jets with a moderate 
Lorentz factor F = 5, corresponding to relativistic though not ultra-relativistic jets, and fol- 
lowed the space-time evolution of the plasma-jet configuration. These simulations are well 
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Fig. 12 Plasma density structure in two 3D PIC simulations with 
different a = and t7 = 0.1 for two colliding pair plasmas with 
F = 15 (data taken from Spitkovsky 20051. Top: In the absence 
of initial magnetic fields the Weibel filamentation instability gen- 
erates stretched current and density filaments and the correspond- 
ing transverse (toroidal) magnetic fields. This causes a comparably 
broad shock transition region ultimately yielding a shock of com- 
pression ratio Ni/Ni fs 7/(7— 1) and a compressed field which is 
essentially restricted to the shock, whose width equals the skin depth 
of ra UfiXg. Downstream the plasma is thermal but particles of gy- 
roradius larger than the shock width can escape upstream because 
they move at speed c, while the shock moves upstream with speed 
c/3 only. Bottom: In the magnetised case, the upstream flow runs 
into the shock front, is retarded and causes a narrow transition by 
the diamagnetic current of stopped and gyrating particles as well as 
a stationary precursor wave which is radiated from the shock. The 
ramp cuiTent increases the magnetic field over a length of < \OXe 
or a few gyro-radii in the compressed field, leaving it enhanced and 
turbulent over some downstream distance. Hence, non-magnetic and 
magnetic shocks obey a completely different physics and have a 
completely different structure; they are controlled by different in- 
stabilities. Since, however, the magnetic field is generated in a very 
short time, any ultra-relativistic shock will readily become magnetic 
after a short time. 



suited to investigate the mixing of fast jet material with the environment by effects which 
are caused by the presence of the lateral jet boundaries. What concerns the formation of the 
relativistic shock which is generated at the head of the jet so the main result was the genera- 
tion of the Weibel-filamentation instability which was observed in all cases and found to be 
strongest for an initial ct = 0. 

Effects of Magnetisation. The cr-dependence of the relativistic shock properties have been 
investigated more closely by ,Spitkovsky^ ( ,2005 j in a number of simulations of pair plas- 
mas in head-on collision of two counter- straming quasi-neutral plasmas. This arrangement 
is perfectly suited for the investigation of relativistic shock formation and monitoring the 
evolution of shock structure that is independent of the effects of the distant lateral plasma 
boundaries. 

The most important finding that could be extracted from these simulations was that the 
shock structure indeed sensitively depends on ct. At low magnetisations ct < 10^^ the shock 
turns out to be indeed mediated by the Weibel instability. However, at larger magnetisation, 
i.e. in the presence of external fields or also after the Weibel fields have run nonlinear, the 
shock couples with other waves (such as Bell modes) and settles in high magnetic field 
saturation. The shock is then completely controlled by magnetic reflection of particles and 
substantially less by the Weibel instability. Apparently the transition is gradual with increas- 
ing ct; the newly generated magnetic field gradually takes over and starts dominating shock 
dynamics via the Lorentz force that it exerts on the particles. 

What concerns the downstream particle distributions, [Spitkovsky] ( |2005| l argues that 
these distributions are always thermal, apparently for all ct, though this claim is not yet 
proved by the simulations themselves and might in its stringency possibly hold only for pair 
shocks. In fact, the simulations by |Spitkovsky| l (2005| l show indeed no direct sign of any 
shock-particle acceleration. But we may recall that the ID-simulations by |Amato and Arons] 
( [2006| ) as well yielded a thermalised downstream distribution, but they also indicated the 
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Fig. 13 Left: The paths of four selected particles (in colour) out of the nonthermal tail of the downstream 
particle distribution drawn on top of the magnetic field as function of time. The particles manage to pass the 
shock ramp being reflected and accelerated. The gain in energy is shown in the right part of this panel for all 
four particles. Inspection shows that the main energy gain is accumulated in resting at the shock ramp and 
moving back and forth around the ramp on small oscillations only. Right: The particle spectrum downstream 
in a lOOAf wide slice at downstream distance 500/l(,. Maxwellian fits are shown, and a power law fit to the 
flat region is indicated, (data taken from |Spitkovsky|2()08b^ . 

presence of a small number of quite energetic particles in both downstream and upstream 
directions. After a sufficiently long time those among the latter having upstream directed 
velocities might in nature outrun the slightly slower shock and reach the shock front, possi- 
bly even escaping to upstream to be further accelerated. We might note that any positively 
charged downstream particle that makes it up to the shock front experiences the upstream 
accelerated shock potential which may kick it out into the upstream shock vicinity. Never- 
theless, these simulations of relativistic shocks do not unambiguously support the view that 
shock acceleration is acting and is indeed the canonical acceleration mechanism in the Uni- 
verse. |Spitovs^([2005]l argued that there are probably no shocks with low-magnetic field 
with dimensions > ID in nature. 

Long-term Behaviour. In order to check the long-term behaviour of shocks, |Spitkovsky| 
(p008b ) also performed 2D PIC simulations in unmagnetized pair plasma and found that 
the particle distribution downstream of such shocks consists of two components: a relativis- 
tic Maxwellian, with a characteristic temperature set by the upstream kinetic energy of the 
flow, and a high-energy tail that extends to energies > 100 times that of the thermal peak. 

This high-energy tail is best fitted by a negative power law in energy of index 2.4 ± 0. 1 , 
modified by an exponential cutoff owing to the finite size of the simulation box ( [Bykov| 
[and Uvarov||1999| l. The cutoff moves to higher energies with simulation time, causing an 
increase of the power-law range. The number of particles in the tail is approximately esti- 
mated as Nfaii ~ 0.0^2. However, these particles carry a substantial fraction of w 10% of 
the downstream kinetic energy. The idea is that particles from the hot thermal downstream 
background above a certain energy with velocity directed upstream are in the long term in- 
deed fast enough to outrun the comparably slow shock which moves with nominal speed 
As/i ~ c/3 upstream. Once in the upstream region, these particles are believed to be back 
-scattered and to undergo Fermi-type acceleration. 
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Fig. 14 Left: Simulated spectrum of the 2D-magnetic power P(k) and of the evolution of the fj parameter 
taken in the downstream interval — lOOOAc < Ax < — 200Ae at three progressing times, with k = 27l/X wave 
number, A .x downstream distance from the shock ramp, and A,, = caig the electron skin depth. The increase 
of the integrated magnetic energy fraction (7(t) with time T is shown as big dots. The diagonal corresponds 
approximately to the demarcation line A = fv between small-scale (magnetised, wavelength < rms-gyroradius 
rc ~ v^2P(i:)Ae) and large-scale (non-magnetised) fluctuations. The maximum of the growing power is for 
all times above this line indicating maximum growth of downstream fluctuations at small wave numbers-long 
wavelengths (i.e. at larger scales). Artificially cooling the system by replacing pailicles above F = 80 with 
thermal particles stops the wave growth as indicated by the pink ('cooling') curve. Right top: Normalised 
electromagnetic energy density for two simulation times Ti (O^ = 2250, T2 0(, = 1 1925, early (blue colour) 
and late (red colour) times, respectively. Note the increase towards higher fluctuation levels at larger distance 
to both sides of the shock. Right middle: Normalised density N/Nq. With time the shock profile flattens a 
bit towards downstream. Right bottom: Average pailicle momentum taking all particles and only downward 
streaming particles, respectively. Note the substantial increase in upstream energetic particles at late time. The 
upstream flow starts retards substantially closer to the shock ramp than one observes the upstream energetic 
particles. Near and behind the shock the main momentum contribution comes from the energetic particles 
(after jKeshet et al|2009f . These simulations show the presence of an over-abundant population of high-energy 
accelerated particles upstream of the shock. 



Spitkovsky ( 2008b i analysed the trajectories of some of these particles in a similar way 
as |Jaroschek et al| 1 2004| |2005[ ) suggesting that the particles indeed may have escaped up- 
stream and were back-scattered in the magnetic fields generated self-consistently by the 
Weibel instability. If a large enough number of such particles can accumulate with time in 
the tail of the downstream distribution they could in addition enhance the backscattering 
efficiency by exciting Bell modes and their oblique relatives. Subsequently the shock would 
become strongly magnetised, possess a magnetised foreshock and would become mediated 
by the self-generated energetic particle component. 

Because of their small boxes and comparably short simulation times, all these simu- 
lations, even those of |Spitkovsky| l (2008b| l, are incapable of detecting any massive shock 
acceleration of particles, though we already noted some downstream indications for the po- 
tential presence of accelerated particles whose origin still remains unclear; it might probably 
rather to be searched in a process of downstream turbulent acceleration than in a Fermi-like 
mechanism. 
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Recently [Keshet et al| ( |2009[ l extended both the length L of the box and the simulation 
time T up to a combined number (L/Ae)^Ta)e in non-magnetic (initial ct = 0) pair-plasma 
shock simulations with F = 15. They used the same code as [Spitkovsky | p005| but replaced 
the two interacting plasma streams by one stream only and applied a reflecting-wall set-up 
in a box of total scale 6300^ x 1024 = 4 x 10'° cells. The total number of particles was 
increased to w 2 x lO'". 

What interests us here is the long-term evolution of the shock (i.e. its state close to the 
end of the simulation time TCOg > 10"^) as the short term evolution is identical to that de- 
scribed above. The important findings in these numerical experiments can be summarised 
in four main points: (i) The gradual generation of 'larger scale' and 'larger amplitude' 
magnetic field fluctuations surrounding the shock transition - the term 'larger' will be 
defined more precisely below - than in the former short-term simulations. They also ex- 
tend far into upstream, possibly (though not necessarily) suggesting that energetic parti- 
cles manage to pass the shock front from downstream to upstream, where they contribute 
to the generation of magnetic fluctuations. One should, however, note that these authors 
do not follow single particle orbits and thus cannot infer about the origin of the high- 
energy particles, (ii) An increase of the total normalised electromagnetic fluctuation en- 
ergy e^,,, = a + {\^\^) / meNic^ (Fq — 1) occurs both downstream and upstream of the shock 
ramp.(iii) The density profile of the shock is modified by the presence of energetic accel- 
erated particles, and the shock slightly accelerates. This is an effect of the presence of the 
increasingly hotter (more energetic) downstream plasma, (iv) The upstream energetic parti- 
cles increase the average upstream momentum far ahead of the shock at distances ~ 10^^ A;, 
while the bulk momentum of the upstream flow starts decreasing only at upstream shock- 
distances < lO^Ag. 

These results can be read from Figure [14] The left panel shows the evolution of the 
magnetic power spectrum and average a. One finds that the maximum magnetic power 
is generated at 'non-magnetised' scales (larger than the gyro-radius) of the order of A ~ 
3 X which is still very short, yet one order of magnitude larger than the typical Weibel 
wavelengths, and the wavelength does slowly increase with time. Still, the simulations stop 
after a comparably short time and practically no saturation is reached. Also, though the 
wave power increases by somewhat more than two orders of magnitude, its peak in these 
simulations remains about constant at 0,nax ~ 7% at the shock front. What concerns the 
high-energy particles that possibly escape to upstream one should, however, carefully take 
into account that the simulations by |Spitkovs^ ( |2008b[ l indicate that this happens only 
in the unmagnetised case while the magnetisation inhibits such an escape for at least a 
while. Figure [14] shows that there is substantial magnetisation caused. Hence, it cannot be 
concluded that the upstream particles are those which have escaped from downstream. For 
this to know the particle orbits should have been followed, also, one would wish to have 
magnetised simulations for long enough times to see whether then sufficiently hot particles 
can be generated downstream for escape to upstream as well and how long this will take. 
In any case, the simulation times are short compared with any natural time scale for direct 
comparison. 

Similar long-term simulations in unmagnetised electron positron plasma were recently 
performed in 3D by |Nishikawa and et al.|f20 10b ) who injected F = 15 and F = 100 electron- 
positron jets into a long though narrow box (on a 4005 x 131^ grid). Their goal was to ex- 
amine the nonlinear stage of the Weibel instability and its capacity for particle acceleration. 
These authors found cold jet and ambient electrons to become thermalized in the result- 
ing shocks. In these simulations the pair jet interacts with the background pair plasma and 
causes an external bow shock. Large amplitude electromagnetic fields generated behind this 
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Fig. 15 The evolution of the 
downstream particle-energy spec- 
trum as function of the upstream 
magnetic inclination angle 6b„ for 
Tfflc = 9000 and T running from 
1 to 2 X 10^^. Below the critical in- 
clination angle Q « 34° particles 
can escape to the upstream region, 
undergo SDA and return to down- 
stream, where they contribute to 
a nonthermal tail on the particle 
energy spectrum. Above the crit- 
ical inclination, no particles es- 
cape to upstream and no super- 
thermal tail evolves (after [Sironi| 
[and Spitkovsky|2009a) . 



shock lead to turbulent downstream magnetic fields. The Weibel instability is indeed iden- 
tified; it evolves nonlinearly to large amplitudes. The authors argue that the field would be 
strong enough for causing the synchrotron afterglow emission in Gamma Ray Bursts. 



Magnetic Field Obliqueness and Acceleration 

Dependence on Magnetic Inclination. So far either initially unmagnetised or magnetised 
perpendicular pair shocks have been discussed. An attempt to include oblique shock an- 
gles has been undertaken by [Sironi and Spitkovsky| l |2009a| ) who performed 2.5D and 3D 
PIC simulations in order to investigate the dependence on 0b„, the angle between the up- 
stream magnetic field Bi and the shock normal n. A set up with Lorentz factor F = 15 is 
used with cell number 50000. The plasma skin depth chosen to be = 10 cells, each cell 
containing initially 2+2 particles. Emphasis is put on particle escape to upstream and par- 
ticle acceleration. In or der to control the Weibel instability, the initial magnetisation is set 
to CT = 0.1 > 10-3 (cf. |Spitkovsky|2005 for discussion). In doing so the simulations can 



directly be applied to Pulsar Wind Nebulae (PWNe) and internal shocks in GRBs and AGN 
jets. 



Sub-luminal versus Super-luminal Shocks. Particles can escape upstream only when their 
upstream speed component along the magnetic field is larger than the sub-luminal speed of 
the shock. Otherwise the shock is super-luminal with respect to the particles, and they will 
be swept downstream. In the reflecting-wall frame the critical angle between sub-luminal 
and super-luminal shocks is tanS^n.en/ = Ftan(cos-' with jSj^ the shock speed in the 
upstream frame; this yields dBn.crit ~ 34° in the wall frame d Sironi and Spitkovsky|2009a i. 
Such quasi-parallel relativistic shocks thus will behave similarly to unmagnetised shocks, 
possessing foreshocks and are subject to mediation by the Weibel filamentation instability 
which is the counterpart of the resonant magnetic-pulsation instability which is excited in 
the foreshocks of quasi-parallel non-relativistic shocks. 

Figure [15] shows the simulated downstream particle-energy spectra as function of F 
and in dependence on the magnetic inclination angle 6 taken at identical simulation time 
TCOe = 9000. The simulated energy spectra confirm indeed that only sub-luminal shocks let 
downstream particles return to upstream and become accelerated. Such particles, in addition. 
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contribute to the generation of upstream waves which, in the relativistic domain, parallels 
what is known from non-relativistic oblique shocks (cf. [Balogh and Treumann||20l"T} . In 
contrast, 'super-luminal' shock^^do not permit particles to escape upstream and therefore 
lack any upstream waves. One may thus conclude that the presence of an oblique upstream 
magnetic field at angles below the critical angle suppresses the escape of downstream parti- 
cles to upstream and thus also suppresses the self-consistent generation of upstream waves 
and further shock acceleration of the escaping particles to higher energies either by the dif- 
fusive shock (DSA) or shock drift (SDA) acceleration mechanisms. This conclusion holds 
for the self-consistent acceleration and upstream injection of particles from the background 
pair plasma and should be valid until the downstream turbulence accelerates particles to en- 
ergies such high that their gyroradii become so large that they effectively demagnetise in the 
downstream plasma. 

In the sub-luminal case the downstream particles which are believed to escape upstream 
become accelerated by the SDA mechanism, finally return to downstream and contribute to 
the non-thermal tail of the downstream particle energy spectrum. Like in Figure [13] the 
downstream particle spectrum consists of a relativistic thermal background plus a non- 
thermal tail with power index lying between 2.8 ±0.1<a;<3±0.1 and having a high- 
energy exponential cut-off. The evolution of the particle spectra for different magnetic in- 
clination angles is shown in Figure [15] The self-consistent DSA mechanism with shock 
generated particle injection works only for nearly parallel shocks. 

4.2 Electron-Ion Plasmas 

Electron-positron (i.e. pair) plasmas may not be rare in astrophysics. They are found in var- 
ious objects, particularly in pulsar winds, pulsar wind nebulae and possibly also in Gamma 
Ray Bursts, where they may be created by various processes. The choice to use pair plasmas 
in simulations has not been guided by their abundance but by the advantage the mass ratio 
jU = 1 offers. 

Real astrophysical shocks, in particular external shocks, involve ions and, because of this 
reason, are believed to be much stronger cosmic-ray accelerators than pair shocks which, 
however, accelerate leptons and thus are directly subject to synchrotron emission. The main 
interest in external shocks should really be on treating electron-nucleon plasmas with, pos- 
sibly, an admixture of an additional dilute pair component. Such plasmas have jU 2> 1, in 
case of a pure electron-proton plasma jU = 1836. This mass ratio puts the time scales on 
which electrons and nucleons evolve far apart yet does not allow to treat the ions as either 
immobile or as a simple fluid. 

In the following we briefly review a few contemporary attempts to simulate relativistic 
coUisionless shocks that include nucleons. Basically two different types of simulations have 
been performed: (i) those where the focus is on the generation (and injection) of energetic 
particles by the relativistic shock-production process as well as on the production of mag- 
netic fields; this resembles the pair plasma simulations; (ii) simulations where a high-energy 
cosmic ray component is assumed to be present and the focus is on its further acceleration by 
considering a beam-plasma interaction which produces strong magnetic fields and scatters 
the Cosmic Rays. In addition, cosmic-ray-modified and also radiation-modified shocks have 

To avoid confusion, one sliould keep in mind tliat the term 'super-luminal' does not mean that the shock 
is moving at velocity faster than the speed of light. It is just the geometric point of intersection of the inclined 
upstream magnetic field and the shock front which moves at an appaiently super-luminal speed along the 
shock. 
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been treated. However, in collisionless shock theory on which we focus, these modifications 
are of secondary importance. Though of interest in cases where the cosmic ray component 
is externally given, we may spare any efforts to look into such modifications until a genuine 
understanding of relativistic shock formation, shock structure and shock properties under 
collisionless conditions has been obtained, a state from which we are still far away. 



No Cosmic Rays 



Initially Non-magnetic Shocks. The first successful PIC simulational attempts by Nishikawa 



et al (2003 IW^, [Frederiksen eTall p004| , [Hededal et al ('2004' "lOOF), and A mato and 



Arons (2006) to infer the 3D-shock structure, magnetic field generation and particle accel- 



eration in nonmagnetised relativistic electron-proton plasmas still used a small mass ratio 
of jU = 16. [Frederiksen et al| ( |2004| ) worked with a Lorentz factor F = 3 and performed 
simulations on a small grid 200^ x 800 with 25 particles per cell, a total of 8 x 10^ particles 
and box sizes (lOA,, 1 OA,-, 40 A,). The shock was generated in counter-streaming (along z) 
plasma interaction. One of the plasma streams was taken to have 3 times higher density, and 
the simulations were performed with periodic boundary conditions in (x,y) and open condi- 
tions in z- The simulation time zcOe = 480 (corresponding to TO),- = 120) was long enough to 
allow the relativistic particles to cross the box along z roughly ~ 3 times. As expected, this 
simulations allows for the electron Weibel filamentation mode which grows nonlinearly and 
deflects ions. At later times the ions react with the ion- Weibel mode. Electrons are heated 
by mixing the two streams. The power spectrum of the magnetic field is a power law toward 
smaller scales, and cascades inversely in the transition region from electron to ion Weibel 
modes. 

Larger Lorentz factors F = 15 which are more suitable for ultra-relativistic shocks, and 
different numbers and sizes ofcells (125^x2000) were used by |Hededal et al| ( |2004||2005| ). 
These allow to examine in detail the energy evolution of turbulently heated electrons. It 
turns out that their spectrum is quasi-thermal yet contains a very short power law region 
just before the exponential cut-off produced by escape-losses due the finite length of the 
simulation box. Though this power law range (being shorter than one decade in energy) 
might not be overwhelmingly convincing, the observed substantial heating of the electrons 
in the ion-electron Weibel-mode interaction indicates that relativistic shocks may indeed 
generate high electron temperatures. In nature these temperatures might possibly reach the 
high values needed to explain the observed synchrotron radiation. 

Larger mass ratios between jU = 16 and /i = 1000 have recently been simulated by 
|Spitkovs^ P008a| l in 2.5D dimensions at F = 15. He produced non-magnetised shocks 
in his previously used reflecting-wall PIC code. This time the simulations were run until 
a steady final state was reached. The code uses 10 cells per A;, and thus resolves the ion- 
inertial length scale A,-. The typical stationary shock thickness in this case becomes w 50A/ 
and seems to be practically independent on the mass ratio up to value of jU w 1000. 

As is expected, the dominant mechanism of magnetic field generation can indeed be 
attributed to the ion- Weibel instability. As is also expected, the magnetic energy peaks are 
found right at the shock transition. It comes close to equi-partition with very high final mag- 
netisation ratio CT w 1 locally, but the average magnetic energy remains < 15% of equiparti- 
tion at the shock. The high local value of a is caused by the large overshoot in the magnetic 
field in the shock ramp. 

The downstream region contains isotropised flows and highly structured magnetic is- 
lands. Upstream of the shock the Weibel filaments are non-stationary. Stationarity is reached 
by electron heating which stops the ion- Weibel instability from growing. Electron heating is 
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Fig. 16 Left: Time-evolution of tlie orbits of tlie 80 most energetic ions in a non-magnetised relativistic 
sliock simulation with F = 20. The particles are coming from the upstream flow, are back-scattered in the 
magnetic turbulence in the shock transition, staying within the distance of an ion inertial length Xj 50Aj, 
and are locally fast enough to traverse the shock a short distance back upstream. They do not make it far 
into upstream but oscillate around the shock thereby picking up energy as shown in the inset for one particle. 
One should note that all particles during their shock-drift phases remain within the ion-inertial length shock- 
transition region of Ax ~ 50A(, cross-shock extension, as indicated by the shading thus moving together with 
the shock upstream at the shock velocity. The colour code along the orbits of the particles indicates the 
increase in particle energy. Once the particles have gained sufficient energy to overcome the reflecting shock 
potential and the scattering upstream turbulence, they escape with the flow into downstream direction where 
they produce the high energy tail of the particle spectrum. Right: The increase in the proper ion Lorentz 
factor along the particle orbits. The energies along the orbits form a cloud showing the average increase of 
the Lorentz factor for all 80 ions up to saturation at about 150 (af'ter |Martins et al|2009b^ . 



large with average downstream-electron energy ratio e^. ~ 0.5 as needed in GRB afterglows, 
but the electrons are mainly thermal. Hence, it seems that electron-proton shocks behave 
rather similarly to relativistic pair shocks. 

Similar results have been obtained by [Martins et al| ( [2009b[ l for a relativistic electron- 
ion stationary shock in a PIC simulation with initial F w 20 but small fi = 32 and with a 
resolution of 10 cells per A^,. This high resolution was chosen in order to follow the single 
particle dynamics. As before, the magnetisation evolved up to ~ 20% of equipartition near 
the shock, somewhat less than in the above described simulations by |Spitkovsky| l |2008a| ) 
which is probably due to the smaller mass ratio. Interestingly, both electron and ion energy 
distributions evolved into high energy tails up to proper Lorentz factors y > 100. These 
authors picked out the 80 most energetic particles in their simulation and followed their 
orbits. Those energetic ions who became most strongly accelerated started hanging around at 
the shock front for several 100 plasma times in which they picked up energy until being push 
to their energy limits and being able to overcome the reflecting shock potential and make it 
to downstream. This process is very similar to shock-foot formation in quasi-perpendicular 
non-relativistic shocks with the main difference that in the relativistic case the magnetic 
shock field is self-consistently generated by the Weibel-mode instability. This is shown in 
Figure [16] It is most interesting that their is no indication for the ions to run the shock out. 
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Instead, all the high energy ions seem to originate in the upstream region and do for quite 
long times not escape from the shock transition staying within an ion skin-depth distance of 
the shock ramp. 

Magnetised Shocks. Since astrophysical shocks in PWNe, GRBs and AGNs are magnetised 
(as is obvious from observation of synchrotron radiation), it is important to investigate the 
effect of magnetic fields on shock formation. In non-relativistic shocks magnetic fields are 
the most important ingredient as well, because no collisionless non-relativistic shocks would 
exist without them (we do not speak of the microscopic electrostatic shocks here). 

The dependence of PWNe termination shocks and internal shocks in GRBs on the pres- 
ence of an ambient magnetic field has recently been subject to 2.5D PIC simulations in 
electron-proton plasmas ( [Sironi and Spitkovsky|2010) i for small mass ratios 1 < /i < 100. It 
turns out that, like in the pair plasma case, relativistic ions (i.e. protons) can return along the 
magnetic field and enter the shock upstream region only when the shock stays 'subluminal'. 
This implies that the shock must be nearly parallel with magnetic inclination angle 0g„ < 30° 
for F > 5 and a > 0.03. Such ions are further accelerated by 'surfing' along the shock. 
Their downstream energy spectrum evolves into a thermal spectrum with an energetic-ion 
power-law tail containing ~ 3% of the downstream ions while consuming ~ 10% of flow 
energy. Electrons are accelerated much less strongly. This is due to their stronger magnetisa- 
tion. Their number is one order of magnitude smaller than that of the energetic ions. These 
findings are not particularly encouraging concerning the ultra-relativistic shock particle (by 
either SDA or DS A) acceleration, for it must be expected that most magnetised collisionless 
shocks are very close to being perpendicular with Qgn ^ 80° . 

These simulations are supported by an independent 3D PIC simulations ( |Hededal and| 
|Nishikawa|2005| l where a relativistic plasma jet was injected into a plasma immersed into 
an ambient magnetic field either parallel or perpendicular to the jet. The direction of the 
ambient field affects the evolution of the electron Weibel instability and the associated shock 
in simulations with (low) mass ratio ^ = 20 when the jet has Lorentz factor F = 5. Like in a 
pure pair plasma (as discussed above), the Weibel two-stream instability grows in the parallel 
and weakly oblique case below the threshold magnetic field strength a ~ 0.03 where it gives 
rise to a transverse magnetic field. 

If the ambient field is perpendicular, the picture in the jet case becomes very complex. 
Behind the jet front (i.e. the shock) inertial separation between electrons and ions causes 
charge-separation fields in which both ambient and jet electrons become strongly accel- 
erated. This kind of electrostatic acceleration, or heating of the electrons, transfers excess 
energy stored in the ions to the electrons and thus provides a way of electron accelera- 
tion which is, however, restricted just to the shock front. In an extended jet-shock front the 
electrons in this field remain magnetised and will in addition perform an E x B drift along 
the shock front, which gives rise to an electron sheet current, whose magnetic field is high 
while being restricted to the shock where it generates the shock overshoot field. This cur- 
rent can become quite strong, probably even strong enough to lead to the excitation of the 
Buneman instability (or to the modified-two stream instability), which both result in the for- 
mation of electron holes in the shock transition region and cause further electron heating 
and non-thermal acceleration. These mechanisms have not yet been sufficiently explored 
for relativistic shock though being well-known in non-relativistic shock dynamics ( |Balogh| 
[and Treumann'201 1). First investigations in this direction have only recently been published 
jbieckmann 2009, Dieckmann et a l,20I0[ [Murphy et al|20I0| . In particular, [Murphy et al| 



2010 1 in weakly relativistic numerical simulations of a shock, found that the shock evolved 



into a quasi-perpendicular shock just because of the presence of this drift current sheet which 
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became unstable and formed magnetic structures, i.e. magnetic flux tubes of fairly homoge- 
neous strong magnetic field. Relativistic electrons gyrating in these field may contribute to 
the generation of synchrotron radiation. 



Cosmic Ray Beams and External Shocks 

Simulations including Cosmic Rays put emphasis on the generation/amplification of mag- 
netic fields either by the Weibel filamentation or Bell-like modes (cf., e.g., |Amato and Arons| 
|2006[ [Amato et al|2008| l, the modification of the shock, increase of the compression factor 
and further acceleration of Cosmic Rays, boosting their energies up as much as possible. 
This kind of approach is less suited for understanding shock formation than for purposes 
of application (mostly to external SNR shocks) where the shock does already exist. The 
cosmic-ray particles are then either generated in the self-consistent and still badly under- 
stood evolution of the shock or - and more probably - originate from external sources, en- 
counter the shock and are accelerated to higher energies, modify the shock and contributinge 
to amplification of the shock magnetic field. In the following we only consider recent simu- 
lations where cosmic ray beam plasma interactions are investigated in view of generation of 
magnetic fields (Weibel and Bell modes). 

Typical approaches of this kind are found in |Niemiec et al|p008[ l who, in view of appli- 
cation to the external shocks found in SNRs, inject a weakly relativistic F = 2 cold cosmic 
ray beam of current Ja- = eNc-cPsh into the plasma which is assumed to drift with the shock 
along the magnetic field So|| into the thermal upstream medium. jSj/, is the normalised rel- 
ativistic shock speed. The shock is assumed as quasi-parallel. The set-up is charge neutral, 
and the current is compensated by the return current 7^, = —eNeVii carried by the back- 
ground electron component. Periodic boundary conditions and various 3D and 2D grids 
have been used. Mass ratios range from = 3 to /i = 500, and the angle of the ambient 
magnetic field is varied from = 40° to 80°. The main result of these simulations is that the 
growth of the magnetic field is slower, strongly affected by the obliqueness of the field, and 
the saturation level substantially smaller than analytically predicted by linear and non-linear 
theories and by the assumed saturation through magnetic particle trapping and scattering (as 
has been discussed in earlier sections), which would just yield (|b^|)/SQ|j w 1. Regarding 
the effect of the magnetic field turbulence on the particle distribution it is conjectured that 
in all simulations alignment of the bulk and cosmic ray flows is obtained. 

Additional simulations have recently been performed by the same authors with cosmic 
ray beam-plasma interactions at much higher Lorentz factors up to F = 300 in 2.5D with 
mass ratio /i = 20 and cosmic-ray flow along the ambient magnetic field ( jNiemiec et al| 
|2010| l. Referring to the linear growth rate of the Weibel or Bell instabilities, the simulation 
time allows for 20 to 30 e-foldings of the growing waves. These simulations essentially con- 
firm the above result of generation of magnetic fields bei either instability while reaching 
much higher field amplifications up to (|b^ |)/5g|| w 5 x 10^ in the highest ultra-relativistic 
case and maximum-growing wavelengths with ?imax/h approximately 10 to 30. It is, how- 
ever, not entirely clear what this large amplification factor actually means if the simulation 
starts with initially zero magnetic field. The amplification is roughly proportional to F, with 
the parallel field being less amplified by about one order of magnitude. The reason for the 
magnetic field amplification is the non-resonant interaction between the cosmic ray beam 
and the plasma, and the Cosmic Rays respond by being scattered in the wave field approxi- 
mately corresponding to Bohm diffusion. 
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4.3 Radiation 

Self-consistent Photon Spectra from Simulations. 

The ultimate link (and proof of astrophysical validity) between the theoretical and simula- 
tion results is the reproduction of observed radiation spectra by radiation spectra inferred 
from simulated shocks. Until recently, such attempts have not been particularly success- 
ful, because the 'synthetic' spectra from simulations could not be measured directly. With 
the available grid size sufficiently short wavelengths, which occur in real space, cannot be 
resolved. Spectra were inferred subsequently from the particle distribution by applying ra- 
diation theory, which is a reasonable approach as long as radiative losses and the shock 
response to the radiated photons can be neglected. 

Recently a different approach has been developed which is based on post-simulation pro- 
cessing codes which have become available. With their help radiation spectra can be calcu- 
lated directly from the particle motion by considering the sub-gyro-orbit 'jitter'-acceleration- 
deceleration sequences of electrons and positrons when following the particle orbits and by 
use of the far field representation of the radiation field ( |Hededal|2005[ [Martins et al|2009a| ). 
A theoretical account for this radiation including simulation results has been given recently 
l [MedvedeV|2009b , Medvedev et al|2010[ l. The new methods allow for resolution up to 1% 
of the cell size and have been applied to both collisionless pair and electron-proton shocks in 
order to determine the shape of the emitted photon spectrum ( |Sironi and Spitkovsky|2009b[ 
[Nishikawa and et al.|20 10b a; Frederiks en et al|2010[ l. 

Sironi and Spitkovsky| (|2009b| exploit ~ 10"* energetic electrons moving in the self- 
consistently produced fields of their 2D PIC shock simulations and were followed over 135 
plasma times. Such spectra must be regarded as local and instantaneous. The photon spectra 
obtained in these simulations obey the expected shape of synchrotron spectra of electrons 
moving in a strong magnetic field. Toward low frequencies all calculated spectra consis- 
tently decay like a ~ co^J power law, in perfect agreement with the expectation of a 2D 
synchrotron spectrum emitted in magnetised media. Qualitatively though not quantitatively 
(yet not surprisingly) this result confirms that the downstream electrons, which have been 
self-consistently accelerated in the downstream shock turbulence, are the sources of syn- 
chrotron radiation that is emitted by relativistic shocks. 

[Frederiksen et al| l |2010[ l, in contrast, in a similar investigation found a number of addi- 
tional results in 2D and 3D PIC simulations of pair and electron-proton plasmas, the most 
prominent one is that the spectral peak increases first as cOpeak °= and then as<x F^. More- 
over, the spectral shape and peak depends on angle. The spectra are anisotropic, indicating 
beaming, which seems more pronounced in electron-proton than in pair plasma. Their spec- 
tra, in contrast to [Sironi and Spitkovsky| ( |2009b[ l are interpreted as non- synchrotron and 
become synchrotron only if the magnetic field is taken as static and not self-consistent. In 
addition, their spectra differ substantially in 2D and 3D simulations. 

These methods, as ingenious as they are and as intriguing and seducing as their results 
may be, should however be taken with substantial care. Not only that the two above men- 
tioned calculations produced different and partly contradicting spectra. The spectra derived 
from the simulations by these methods are by no means real particle spectra but spectra ex- 
tracted from macro-particles which do not account for the emission of each single electron. 
They, instead, refer to a very large coherent clump of (non-physical) particles with identical 
or hidden internal dynamics, and it is not clear what the radiation spectra which are deter- 
mined in the above described way really mean. Their uncritical application to real problems 
must therefore be cautioned. 
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Radiation Mediated Shocks. 

Another broad research field in relativistic shocks is centred on the question to what extent 
shocks are modified by their own self-consistently produced radiation. Radiation may in- 
deed modify shocks, but the effect is different from that caused by the production of Cosmic 
Rays. Radiation results in energy loss and, if reabsorbed, causes redistribution of energy in 
the shock environment. This problem becomes urgent in hot relativistic shocks which emit 
very intense radiation. Such shocks have so far been excluded from our discussion. Here, we 
just note that |Budnik et al| ( (2010[ ) have recently investigated some of the consequences that 
intense self-generated radiation has on the relativistic shock structure for Lorentz factors in 
the interval 6 < F < 30 under the assumption of coupling of electrons and protons through 
plasma processes. These authors provide the relevant equations and boundary conditions 
for later application. Shocks of this kind become radiation-mediated. Simulation of the re- 
lated problems would involve the self-generated radiation field in the equations of motion of 
the simulated macro-particles, which also includes pair-production in the intense radiation 
fields. So far the radiation produced in simulations is, however, undetectable. In order to 
observe its effect its intensity would need to be up-scaled to the expected local strengths. 
This is still far from any practical application. 

Relevance for Modelling Gamma Ray Bursts. 

Following the original suggestion by [Remington et al| ( [1999| l, [Medvedev and Spitkovsky] 
l |2009aj ) have addressed the important question, whether such 3D PIC simulations (and also 
related laboratory laser plasma experiments like those performed by |Kuramitsu et al|20Il| ) 
are just qualitative and, quantitatively, will or will not in the near future become comparable 
with astrophysical observations, as for instance the radiation from GRBs. In order to answer 
this question, one needs to compare the radiative cooling time of relativistic electrons in the 
self-consistently generated shock magnetic field and the microscopic dynamical time of the 
evolution of collisionless relativistic shocks, i.e. the inverse proton plasma time Tcooi(Oi- In 
current simulations in 3D, the result of such a comparison is that, for Tcooi(i>i ^ 10^, electron 
cooling in the vicinity of the shock becomes efficient such that the electrons contribute 
to a substantial loss of kinetic energy in the form of emitted radiation. Current 3D PIC 
simulations can, indeed, resolve this type of relativistic shocks, and thus permit the detailed 
study of their micro-scale structure. 

With modem methods (and the caution noted above) it might become possible to infer 
the radiation spectrum of the electrons, and even to determining their spectral slopes and 
peak or integrated powers. The radiation is, however, restricted to that emitted from the 
close vicinity of the shock corresponding to the motion of the electrons in the simulation. 
Such spectra can, with caution, be compared with observed radiation spectra to provide 
information about the particle population, energy, and the type of shock responsible for the 
generation of the radiating particle distribution. Shocks of this kind are found as internal 
radiative shocks in baryon dominated GRBs. On the other hand, radiative cooling limits any 
Fermi acceleration (if present) of the electrons, and this implies that the electron spectrum 
is thermal, i.e. lacking a high energy tail, and should peak in the multi-MeV range. In the 
opposite case when the radiation efficiency of the electrons is less than that predicted by the 
simulations, the electron distributions may possess energetic tails. External shocks, on the 
other hand, like those believed to be responsible for the extremely early afterglow phase, 
will probably be cooled only when they are still highly relativistic. 
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5 Conclusions 

Despite intense and increasing research over the last decade, which is reflected in the vastly 
growing number of recent publications on this subject, understanding of relativistic colli- 
sionless shocks is still far from completion or from a situation, where derivation from first 
principles becomes feasible. The reason for this lack of understanding lies in the practical 
non-availability of collisionless relativistic shocks, neither on Earth nor in the near-Earth en- 
vironment though this might change ( [Remington et al|1999| l when laboratory experiments 
will manage to generate highly relativistic shocks in laser fusion plasmas (cf. the comments 
and references given in |Medvedev and Spitkovsky|2009aj i. In particular the three questions 
of interest in astrophysics, (i) the generation of strong large-scale magnetic fields in shock 
interaction, (ii) the related shock acceleration of Cosmic Rays to very high energies, and 
(iii) the self-consistent generation of radiation still lack any sufficient understanding and 
mapping to observations. Yet progress cannot be denied and has been accounted for in this 
review (and also the in-depth discussion given by |Medvedev|2009a| i. 

In the following we very briefly summarise what, as seen by us, is the current state of 
the art with respect to a final answer to these questions. 

Magnetic Field Generation 

Generation of magnetic fields in relativistic shocks is believed necessary and unavoidable 
in increasing the shock compression ratio, as required by the DSA mechanism (in case it 
will be confirmed that this is the only, unavoidable, or the dominant mechanism of particle 
acceleration), as well as for the generation of synchrotron radiation. 

The main processes of magnetic-field generation in relativistic shocks that have been put 
forward during the last two decades are based either on the Weibel filamentation instability, 
which has been confirmed by the numerical simulations in pair and electron-ion plasmas, or 
the generation of Bell-like modes which has also been confirmed in numerical simulations. 

The Weibel instability works best in initially non-magnetised (or very weakly magne- 
tised) low-temperature relativistic flows or in flows about parallel to the ambient magnetic 
field. Saturation levels are, however, relatively low, and are restricted to distances of few 
(ion-nucleon) skin depths, not far downstream of the shock. However, the Weibel modes re- 
magnetise the shock transition such that the shock becomes about magnetically perpendic- 
ular. In fact, being a transverse mode, the Weibel filamentation instability generates quasi- 
toroidal transverse magnetic fields by splitting the original homogeneous beam up it into 
beamlets or beam filaments on the fastest growing transverse scale. This implies a highly 
structure magnetic field along the shock surface different from the model of a homogeneous 
quasi-perpendicular shock plane. Since the filamentation instability, however, provides an 
entire spectrum of scales, the structure of the magnetic field along the shock surface be- 
comes complex. 

The finding in numerical simulations that the downstream region carries magnetic fields 
far more distant from the shock transition than the Weibel mechanism predicts is not well 
understood yet. These fields form vortices (as first identified in most recent simulations of 
[Murphy et al|2010l l, which might suggest that in 3D the Weibel fields in the shock transition 
evolve nonlinearly into current loops and (possibly reconnected) field vortice^^ which are 

When this happens, it would, on the other hand, put into question the common assumption that such 
field vortices stabilise the shock transition by particle trapping. At the contrary, one expects that the shock 
transition becomes highly structured and complex in the tangential direction (as suggested by [Murphy et alj 
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advected far downstream by the shocked turbulent flow and thereby possibly fill an extended 
downstream region with small-scale magnetic fields. 

Bell-like modes require the presence of a cosmic-ray population upstream of quasi- 
parallel shocks. This population is either generated self-consistently by the shock (i.e. raising 
the injection problem) or it is imposed from the outside. The Bell-like fields occupy the 
shock-foreshock, causing the foreshock to become turbulent on short scales yet offering the 
possibility of diffusive interaction which drives a secondary wave population possessing a 
longer- wavelength spectrum of scales up to several cosmic-ray gyro-radii. 

The turbulence caused by these waves retards the upstream flow though possibly only 
slightly. Since the turbulent waves have both a very low frequency and a very low upstream 
phase speed, the relativistic upstream flow must advect them toward the shock transition, 
where they accumulate. Depending on how much time this takes, the waves with longer 
wavelength will grow to non-linear amplitudes and contribute to a nonstationarity of the 
shock, thereby forcing the shock to reform on irregular time-scales. This reformation will 
be similar to what is known from upstream magnetic pulsations in non-relativistic quasi- 
parallel shocks ( jBalogh and Treumann|20ir) l, a process which has not been given any atten- 
tion so far in relativistic shock theory! 

Once this happens, the shock will continually radiate Bell-like (and also Weibel) mag- 
netic fields to the downstream region. Being frozen into the flow like quasi-stationary struc- 
tures, magnetic vortices and flux tubes of different scales, these fields become advected 
downstream by the flow with very little turbulent attenuation only thereby contributing to 
and extended region of downstream turbulence. They will gradually occupy almost the en- 
tire downstream domain up to distances very far away from the narrow shock ramp and 
shock transition. Since the plasma on the scales of these structure is manifestly collision- 
less, they indeed experience only very little attenuation. This process might provide the main 
mechanism of generating highly turbulent magnetic fields in the shocked medium far behind 
relativistic collisionless shocks. 



Particle Acceleration and Radiation 



There can be no doubt that relativistic and, in particular, ultra-relativistic shocks accelerate 
charged particles to high energies. In the present review we have not dealt with particle 
acceleration. There is a wide literature on this subject, mainly favouring the DSA mechanism 
the consultation of which the reader is referred to. However, it remains unclear whether the 
DSA mechanism is involved in the acceleration process at ultra-relativistic shocks. Available 
simulations on relativistic shock formation in the absence of external Cosmic Rays provide 
little or at least no unambiguous indication for the DSA mechanism and would restrict it just 
to the case of highly magnetised nearly parallel shocks. 

The problem is buried in the following: in very fast shocks it is difficult for upstream 
directed particles of the thermal downstream particle distribution to outrun the shock. These 
particles must originate in the downstream thermal tail and must have been pre-accelerated 
by some turbulent downstream mechanism to Lorentz factors substantially larger than the 
shock's proper Lorentz factor such that their upstream speed exceeds the shock velocityp^ 



|2010^ , with saturation and stabilisation of the Weibel fields provided by non-linear effects quite different from 
trapping, an effect whose consequences have not been investigated yet. 

''^ It is sometimes argued that in ultra-relativistic shocks also the downstream particles move at velocity of 
light c which is independent of any coordinate frame and, thus, the particles would also move at light velocity 
with respect to the shock system. This, however, holds only for photons. The velocity of any massive particle 
must be transformed according to the rules of relativity, and so, relative to an ultra-relativistic shock, they 
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In addition, for magnetised shocks this outrunning will be made possible only if the shock 
is nearly parallel with magnetic inclination angle against the shock normal, i.e. Ogn < 30°. 
This excludes a large class of relativistic and ultra-relativistic shocks, in particular external 
shocks, which have always been believed to be close to perpendicular. 

If a particle makes it finally up to the shock, simulations seem to indicate that it will 
oscillate for a while around the shock front, experiencing acceleration until it ultimately es- 
capes downstream and contributes to a power law tail. This kind of acceleration resembles 
SDA where the particles during the phases when they enter into upstream close to the shock 
experience the upstream convection electric field and become accelerated in a way that is 
similar to pick-up ion acceleration in the shock-foot region over the part of their upstream 
gyro-orbit exposed to the upstream convection electric field E = — jS x B [ . This may happen 
many times for, when the particles enter the shock from upstream, the strong shock potential 
reflects the particles back upstream and the shock ramp sweeps them along until they have 
picked up sufficient energy to overcome the shock potential and leave from the shock to 
downstream. For this process it is even not at all needed that downstream particles reach the 
shock front. In Figure[T6j for instance, all ultimately accelerated ions have come solely from 
the upstream flow while becoming accelerated to large Lorentz factors F in this process and, 
in principle, the measured final Lorentz factors ~ 150 of the particles in this simulation 
can be taken as mapping the cross-shock potential Ush ^ '«/C^r^„/e (w 150 GV in this par- 
ticular case; accounting for the reduced mass used in the simulations this amounts to w 4.5 
GV in a proton plasma. The simulations use a F = 20 which implies protons of energy ~ 20 
GeV, thus the cross-shock potential consumes roughly 25% of the kinetic energy of the 
upstream flow) which the accelerated particles must overcome. In non-relativistic shocks, 
these particles would be considered as shock-reflected ions, while in the ultra-relativistic 
case with high shock velocities jS.,/, they cannot escape far upstream from the shock. The 
shock sweeps them in front of it upstream over a distance until they can escape downstream 
from the shock. Possibly particle acceleration at ultra-relativistic shocks is just this kind of 
reflection-surfing SDA acceleration. 

It remains unclear, however, how in electron-nucleon plasmas the energy is transferred 
from shock accelerated nucleons to electrons so that the ultimately observed photon-radiation 
spectra are generated. The mechanism of electron acceleration is still badly understood. Re- 
cently, referring to the process of reconnection and generation of strong electric fields known 
from near-Earth space physics, it has been suggested in view of the GRB afterglow that 
electron acceleration could be provided by a manifestly non-Fermi process that is based on 



acceleration in strong local magnetic-field aligned electric potential differences ('Medv edev| 
[and Spitkovsky 2009b ) . Indeed, such local electric potentials being built up, might be the 
ultimate mechanism of electron acceleration (cf., |Treumann an d Terasawa 2001). This pro- 
cess is know from magnetic reconnection (cf., e.g., Treumann et al|2009j 2010) and from the 
strong magnetic field-aligned electric fields in the aurora which most probably are caused by 
shear flows in local turbulence. In the case of the Weibel instability such fields might arise 
from merging of current filaments, while in the shock front they will evolve naturally as a 
consequence of the cross-shock potential and the resulting strong electron cross-field drift 
current in the shock front, which should ultimately disrupt into electron phase space holes 
by the action of the Buneman instability. The accelerated electrons experience a jitter motion 
and cause radiation. We referred to radiation spectra of this kind which were calculated from 
following the jitter motion of electrons and yielded satisfactory results on the shape of the 



will have difficulties to catch up with the shock unless their F factor becomes very high, i.e. unless they have 
already been accelerated to very high energies in the region downstream of the shock. 
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emitted photon spectra. Future refinements of such calculations may indeed reproduce the 
observed spectra, though one should keep in mind the caveat that the simulated particles are 
just macro-particles which do not exhibit the individual motion of the accelerated electrons 
in the real natural objects. The calculated radiation spectra are kind of mock spectra the 
similarity of which to real spectra may be just incidental. Moreover, as long as the electron 
acceleration mechanism remains obscure one may not expect any fast progress. 

In the light of these remarks (cf . also |Medvedev|2009a| i we acknowledge without any 
reservation that the most recent numerical full particle simulations have substantially ad- 
vanced our knowledge about relativistic coUisionless shocks. Undoubtedly, they have con- 
tributed to a much deeper understanding of shock formation, shock structure, and the related 
acceleration and radiation processes. The technical improvements expected in the near fu- 
ture concern a realistic mass ratio /J. = 1840 in 3D simulations of electron-proton plasmas, 
substantial extensions of simulation boxes and simulation times, and increases of cell num- 
bers, particle numbers and resolution. 
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